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1.  Introduction 


The  evolution  of  material  mierostrueture  during  proeessing  plays  a  eommanding  role  in  the 
physieal  properties  and  meehanisms  in  materials.  Henee,  establishing  the  proeessing-strueture 
and  strueture-property  relationships  can  help  accelerate  the  design  and  development  of  new 
metals.  Without  such  relationships,  achieving  meaningful  advances  becomes  a  highly  empirical, 
hit-and-miss  undertaking.  Many  materials  already  have  a  significant  foundational  understanding 
from  decades  of  research,  but  emerging  material  classes  such  as  nanocrystalline  alloys  are  still  not 
as  well  understood.  Experimental  studies  are  time-consuming  and  are  therefore  necessarily 
limited  to  a  select  number  of  compositions.  Computational  modeling  has  the  potential  to  be 
quickly  applied  to  many  alloy  systems,  but  the  results  should  always  be  verified  by  experimental 
data.  It  is  by  marrying  these  2  approaches  that  they  can  inform,  and  accelerate,  one  another.  That 
is,  computational  results  can  be  tuned  by  experimentally-derived  parameters  and  experimental 
decisions  can  be  informed  by  predicted  results. 

Nanocrystalline  alloys  present  an  important  subset  of  metallic  materials.^  Because  of  their  small 
grain  size  (less  than  100  nm),  they  tend  to  possess  extremely  high  strengths,^^^  which  are 
associated  with  grain  boundary  (GB)  strengthening,  i.e.,  the  Hall-Petch  effect.^’^°  Hence,  these 
alloys  can  possibly  enable  new  lightweight  components,  high-strength  conductors,  and  many 
other  applications  that  are  not  achievable  using  conventional  coarse-grained  alloys.  A 
fundamental  limitation  to  their  use  comes  from  their  inherent  thermal  instability,  which  has  been 
extensively  explored  using  both  experiments and  computational  approaches. The  small 
grain  size  produces  an  extremely  large  driving  force  for  grain  growth.  If  not  kept  in  check,  grain 
growth  can  occur  at  modest  temperatures,  even  at  room  temperature  in  pure  materials,  such  as  in 
copper  and  palladium.^^  Because  nanocrystalline  alloys  are  often  produced  via  mechanical 
alloying,^"^’^^  this  grain  growth  phenomenon  provides  a  significant  obstacle  to  consolidating  the 
nanocrystalline  powders  by  traditional  powder  metallurgy  techniques  where  high-temperature 
sintering  is  common. 

There  are  several  methods  that  have  been  devised  to  circumvent  this  temperature  sensitivity  and 
to  stabilize  the  grain  size.  Grain  size  stabilization  commonly  involves  adding  small  quantities  of 
an  insoluble  element  (i.e.,  solute).  Nonequilibrium  processing  (e.g.,  mechanical  alloying,  rapid 
solidification)  is  often  used  to  force  the  solute  into  solution  and,  upon  heating,  the  solute 
segregates  to  interfaces,  such  as  GBs.^^’^^  Alternatively,  the  solute  remains  in  solution  or 
precipitates  out  as  a  second  phase.  These  2  basic  mechanisms  for  stabilization  are  known  as 
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thermodynamic  stabilization  and  kinetic  stabilization,  respectively.^^  Thermodynamically,  the 
solute  is  expected  to  segregate  to  GBs  and  reduce  the  GB  energy  to  minimize  or  eliminate  the 
driving  force  for  growth. Kinetically,  the  solute  acts  to  hinder  GB  mobility  by  diffusion-related 
means,  such  as  solute  drag,"^°  or  by  pinning  boundaries  with  a  fine  dispersion  of  precipitates."^^ 

The  operating  mechanism  of  stabilization  can  be  difficult  to  delineate.  Highly  insoluble  elements 
can  segregate  into  interfaces,  demix  into  discrete  particles,  or  form  intermetallic  compounds  with 
the  solvent.  The  outcome  depends  on  many  different  factors,  such  as  elastic  size  mismatch 
between  the  solvent  and  solute  and  their  enthalpy  of  mixing.  To  achieve  thermodynamic 
stabilization,  solute  precipitation  must  be  kinetically  hindered.^^  Predicting  and  accounting  for 
routes  to  grain  size  stabilization,  then,  is  critical  to  deepening  the  understanding  of  grain  size 
evolution  in  nanocrystalline  alloys. 

This  work  first  focuses  on  the  iron-zirconium  (Fe-Zr)  nanocrystalline  alloys.  In  the 
nanocrystalline  Fe-Zr  system,  experimental  studies^"^’^^  have  found  that  a  nanoscale  grain  size  in 
Fe  could  be  maintained  at  900  °C  by  adding  4  atomic  percent  (at%,  or  just  %  throughout  this 
report)  Zr.  Without  the  Zr  addition,  the  grain  size  in  nanocrystalline  Fe  starts  to  coarsen  near 
400  °C.  Zr  was  chosen  as  a  solute  because  of  its  large  atomic  size  difference  with  Fe  (28%), 
which  drives  segregation.  There  is  also  a  chemical  contribution  through  the  heat  of  mixing 
(-25  kJ  mol~^).  This  negative  heat  of  mixing  indicates  that  it  is  more  favorable  for  Fe  and  Zr  to 
form  bonds  between  one  another  than  to  form  bonds  between  like  atoms.  Hence,  intermetallic 
formation  is  a  likely  outcome  at  high  temperatures.  Intermetallic  formation  is  also  a  likely 
outcome  in  regions  with  high  concentrations  of  Zr  (e.g.,  interfaces).  The  competition  between 
atomic  segregation  and  second  phase  formation  suggests  that  the  Fe-Zr  system  must  be 
considered  from  both  thermodynamic  and  kinetic  perspectives.  In  this  manner,  the  Fe-Zr 
nanocrystalline  system  is  similar  to  many  nanocrystalline  systems.  As  such,  it  is  considered  to  be 
a  representative  nanocrystalline  system  for  the  present  study. 

The  objective  of  the  present  work  is  to  formulate  a  predictive  thermokinetic  model  for  grain 
growth.  This  objective  is  accomplished  by  extending  the  combined,  thermokinetic  model  for  grain 
growth^^’"^^’"^^  to  allow  for  temperature  dependence  of  parameters,  thereby  extending  this  model’s 
predictive  capability  over  a  wide  range  of  temperatures.  The  extended  model  is  used  to  predict 
grain  growth  at  a  variety  of  temperatures  and  solute  concentrations  in  the  Fe-Zr  system.  The 
predictive  modeling  is  then  compared  to  X-ray  diffraction  data  collected  in  situ  during  annealing 
of  nanocrystalline  Fe-Zr  alloys  with  Zr  (atomic)  concentrations  of  0.5%,  1%,  4%,  and  10%.  The 
present  work  shows  the  benefit  of  adding  a  temperature  dependence  to  the  thermokinetic  model  to 
predict  grain  growth  for  this  system. 
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2.  Thermokinetic  Model  for  Grain  Growth 


The  parabolic  kinetics  of  grain  growth  are  related  to  the  presence  of  solute  segregation  to  GBs 
and  the  drag  force  that  results.  In  this  scenario,  the  GB  velocity,  V (=  dD/dt),  is  controlled  by  the 
motion  of  single  atoms  through 


1/  =  ^  (l  -  exp 


where  v  is  the  jump  frequency,  Agbs  is  the  spacing  of  atoms  in  the  GB,  R  is  the  gas  constant,  T  is 
the  temperature,  and  Qd  is  the  activation  energy  for  an  atom  to  transfer  from  1  grain  to  another 
under  a  driving  force  AP.  We  can  assume  that  AP  <<  RT  for  most  cases  and  that  AP  ^  'jh/D 
for  curvature-induced  grain  growth  processes,  where  'jb  is  related  to  the  GB  energy.  Additionally, 
the  pre-exponential  factor  v\gbs/RT  is  related  to  the  interface  mobility  Mq,  leading  to  the 


expression 


dD  7fe  f  Qd 

-  =  M,-exp  -  — 


To  calculate  Mq,  Vooijs  et  al."^"^  showed  that  Mq  =  v\gbs/RT  =  k5/hR  with  5  =3.0  A,  and  h 
and  k  being  the  Planck’s  and  Boltzmann’s  constants,  respectively.  This  gives  a  mobility  of 
0.7518  molm  s“^  that  is  only  dependent  on  the  length  scale  parameter  relating  to  interface 
width,  5.  In  Eq.  2,  both  Qd  and  7^  are  functions  of  the  temperature  T  and  the  solute  concentration 


xo,  giving 


n^"fb{T,Xo)  f  Qd{T,Xo) 

^  =  ^exp - — 


It  can  be  shown  that  Qd  and  %  are  interrelated.  For  instance,  as  solute  segregates  to  the  boundary, 
7b  decreases,  and  there  is  an  associated  increase  in  Qd  because  of  solute  drag  effect.  However,  it  is 
also  important  to  realize  that  there  are  limits  on  Qd  and  7^  related  to  reaching  a  metastable 
equilibrium  when  a  certain  grain  size  D*  is  reached  (when  lb  =  0). 

2.1  Thermodynamics  of  Grain  Growth 

The  thermodynamics  of  grain  growth  is  based  on  the  concept  that  nanocrystalline  materials  can 
be  stabilized  against  grain  growth  by  GB  segregation.  This  concept  is  related  to  the  Gibbs 
adsorption  equation  and  can  be  formally  written  as 


lb  =  7o 


RT\n  (  xq 


STboVA 

D 


+  AHs 
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where  70  is  the  GB  energy  for  the  case  of  a  pure  material,  Vm  is  the  molar  volume  of  the  alloy, 
AHseg  is  the  heat  of  segregation,  Ffeo  is  the  saturated  solute  excess,  and  F;,  is  the  solute  excess  at 
the  boundary.  Note  that  Ff,  is  less  than  F bo  for  the  case  where  the  GB  is  not  saturated  by  solute  and 
that  these  values  are  equal  for  the  GB  saturated  case.  The  term  'jb  is  composed  of  2  terms:  the  first 
reflects  the  GB  energy  in  the  absence  of  solute  (70);  and  the  second  larger  term  modifies  'jb  by 
accounting  for  the  energy  reduction  caused  by  segregation  of  solute. 

Chen  et  al.^^’"^^  have  applied  a  first  order  Taylor  expansion  of  Eq.  4  and  an  equation  of  mass 
conservation  for  nanocrystalline  materials  (Tbo  =  xop{D/Q  +  5))  to  obtain  an  expression  for  70  in 
the  unsaturated  (F^/Ffeo  <  1)  and  saturated  cases  (Ffe/Ffeo  =  1)  that  is  a  function  of  the  grain  size 
D,  i.e.,  ^ 

7fe  =  7i+72|^,  (5) 

where 

7o  ~  Tfeo  [RT  (inxo  —  +  AHsegj  ,  Ffe/F^o 

lo  —  ^bo[RT\nxo  +  AHseg] ,  Tfe/Ffeo 

and 

STboRTVmpS,  Fft/Ffoo  <  1 

SriQRTVm/xQ,  Fb/Ffeo  =  1 

The  limiting  case  in  the  Chen  et  al.  thermokinetic  model  is  that  70  must  be  equal  to  or  greater 
than  zero. 

2.2  Kinetics  of  Grain  Growth 


<  1 
=  1 


(6) 


The  kinetic  aspect  of  grain  growth  is  based  on  the  concept  that  grain  growth  occurs  as  a  result  of 
GB  diffusion.  Diffusion  at  the  boundary  is  affected  by  solute  interactions  with  the  GB  structure. 
Borisov  et  al."^^  derived  the  following  expression  for  GB  energy,  i.e.. 


7b 


In 


^  Dl 
aoA" 


-  Inma  , 


(8) 


where  A  =  1,  a (=  1)  is  for  the  interstitial  mechanism,  oq  is  the  average  interatomic  distance,  5  is 
the  GB  width,  and  rua  is  the  number  of  atomic  layers  forming  the  GB.  The  GB  and  lattice 
diffusion  terms  are  and  Di,  respectively,  and  these  can  be  given  as 

Do  =  r'J.exp  (9) 
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and 


Dl  = 


Qd 

W 


(10) 


where  and  are  the  pre-exponential  diffusion  terms  for  Do  and  Dl-  Chen  et  assume 
a  monolayer  coverage  of  the  GB  (ma  =  1),  r^o  =  l/aa^Na,  and  to  rewrite  %  as  a 

function  of  Qd,  i-e.,  7f,  =  Tbo{QL  —  Qd)-  Hence,  the  driving  force  for  grain  growth  decreases 
with  the  increase  of  Qd-  The  limiting  case  is  when  Ql  =  Qd,  which  results  in  7;,  =  0.  In  a 
similar  form  as  7^,  Qd  can  be  expressed  as  a  function  of  grain  size  D,  i.e., 

Qd  =  Qi—Q2j^,  (11) 


where 

Ql 


and 


Ql  +  RTln  +  {RT  (inxo  -  ^)  +  AH,,g)  ,  r,/r,o  <  1 

Ql  +  RTln  j  -  5^  +  (RTlnxo  +  AHseg) ,  rQTw  =  1 


Q2 


SRTVmpS,  Tfe/Tfeo  <  1 

STboRTVm/xo,  Tfe/Tfeo  =  1 


(13) 


2.3  Thermokinetic  Model 


Chen  et  al.^^’"^^  show  that  by  substituting  Eqs.  5-7  and  Eqs.  11-13  into  Eq.  3,  assuming  that 
exp{—Q2/RTD)  ^  1  -f  Yl'i=i  l-/iK~Q2/RTD)\  and  using  boundary  conditions  of  t  =  0, 

D  =  Dq,  t  =  t,  and  D  =  D{t),  they  arrive  at  an  expression  for  the  evolution  of  grain  size  with 
time,  i.e.. 


1 

2 


(if  +  ^)  (-0  -  -Do)  + 


In 


(  71-D+72  \ 
y  71  ^0-1-72  ) 


=  7iMoexp  {-%)  t, 


(14) 


which  can  be  reduced  to 


1 

2 


{D^ 


D  'i  _i_  r  1 72  Q2  I  1  z'  Q2  ^ 
-t^o)  +  [^  +  RRT  +  2  [rT) 


In 


[  71-D-I-72  \ 

y  71^^0-1-72  J 


(15) 


Interestingly,  the  D'^  —  Dq  term  is  of  a  similar  form  to  the  simple  power  law  (D^  —  Dq  =  2Mo7fet) 
of  Burke  and  Tumbull,"^^  but  the  rest  of  the  equations  terms  stem  from  the  interaction  between  the 
kinetic  and  thermodynamic  aspects  of  grain  growth. 
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Table  1  Thermokinetic  model  variables  and  values 


Variable 

Value 

Units 

R 

8.3145 

Jmol-^  K-3 

h 

6.626  068  X 

2  —1 
m^gs  ^ 

k 

1.380  650  3  X  10-20 

m^gs-^K 

Na 

6.02214129  X  10^3 

mol”^ 

A 

1 

a 

1 

ma 

1 

ao 

8  X  10-^0 

m 

d  (=  maOo) 

8  X  10-^0 

m 

^GBS 

3  X  10-^0 

m 

Do 

5  X  lO-o 

m 

D%lDl 

1 

Mo 

0.7518 

molm  J-^  s-3 

2.4  Model  Implementation 

The  present  model  was  applied  to  the  Fe-Zr  system.  Tables  1  and  2  list  the  constants  and  values 
that  were  used  in  our  initial  calculations  of  the  Fe-Zr  system  to  validate  the  model  with 
experimental  Fe-Zr  results.  Table  1  lists  most  of  the  constants  for  the  model.  Notice  that  many  of 
the  parameters  are  not  sensitive  to  the  type  of  material  system;  several  assumptions  are  made. 
First,  rua  assumes  a  monolayer  of  solute  at  the  GB;  this  may  be  different  in  some  cases.  Second, 
the  initial  grain  size  Dq  is  5  nm.  Third,  the  mobility  term  Mq  is  a  material-independent  parameter 
in  the  present  model. 
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Table  2  Thermokinetic  model  material  parameters 
eonstants  and  values 


Variable 

Value 

Units 

7.1  X  10“® 

m^  mol“^ 

PFe 

7.87  X  10® 

g  m“^  (original  units) 

Zpe 

55.847 

gmol“^  (original  units) 

P  =  PFe/ Zpe 

1.4092  X  10® 

molm“^  (eorreet  units) 

T 

1073 

K 

Xo 

0.04 

70 

0.73 

Jm-2 

AHseg 

92  400 

Jmol“^ 

Ql 

334  800 

Jm-i 

Tfeo 

1.25  X  10“® 

molm“^ 

Table  2  lists  most  of  the  material  parameters  eonstants  for  the  model.  The  upper  part  of  Table  2 
lists  a  few  easily  obtained  material  eonstants:  the  molar  volume  Kn  and  the  material  density  ppe- 
These  eonstants  are  assumed  to  be  equal  to  pure  Fe,  whieh  may  not  be  valid  for  high 
eoneentrations  of  solute.  The  lower  part  of  Table  2  lists  the  state  variables  (temperature  T  and 
eoneentration  xq)  as  well  as  some  of  the  thermokinetie  parameters:  GB  energy  of  pure  Fe  70,  heat 
of  segregation  AHs^g,  activation  energy  in  the  lattice  Ql,  and  the  GB  saturated  solute  excess  F^o. 


The  model  was  implemented  with  the  following  constants  by  setting  Eq.  14  equal  to  zero,  i.e., 


fiD)  =  i 


D  +  ^]  - 
71 


a 


I  72  \ 
0  +  ^) 


+ 


tI  I  72  Q2_  1  (  Q2_t  ^ 


f  71 -P+72 
\  71P0+72 


-  7iMoexp  (-^)  t  =  0 


(16) 


A  numerical  method  for  root  finding  can  now  be  used  to  solve  for  grain  size  D  as  a  function  of 
time  t  in  this  equation.  The  bisection  root  finding  technique  was  employed  herein.  A  brief 
description  of  the  numerical  method  follows.  The  pseudocode  for  the  bisection  root  finding 
technique  is  given  in  Fig.  1,  where  the  function  /  is  equal  to  Eq.  16,  xEower  is  Dq,  and  xUpper  is 
D*  (=  —72/71),  and  the  root  is  equal  to  the  value  of  D  where  /  (D)  =  0.  The  bisection  method  is 
an  iterative  technique  that  assumes  the  root  to  an  equation  /  is  located  within  a  certain  interval 
bounded  by  xEower  and  xUpper.  Hence,  /  (xEower)  and  /  (xUpper)  have  opposite  signs.  If  an 
intermediate  coordinate  is  chosen  (xMid),  the  sign  of  /  (xMid)  narrows  down  the  interval  where 
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the  root  exists  (i.e.,  if  /  (xLower)  and  /  (xMid)  have  opposite  signs,  the  root  lies  in  the  interval 
[xLower,  xMid],  and  if  they  have  the  same  signs,  then  the  root  lies  in  [xMid,  xUpper]).  This 
process  continues  to  subdivide  the  interval  where  the  root  is  until  it  has  been  reduced  to  a 
specified  tolerance,  xTol.  While  other  root  finding  methods  may  be  more  efficient  (e.g..  False 
Position,  Newton-Rhapson,  Secant),  the  present  thermokinetic  model  is  not  expensive  and  the 
bisection  method  converges  to  the  root  within  a  specified  tolerance. 

Input  xLower,  xUpper,  xTol 

yLower  =  f(xLower) 

xMid  =  (xLower  +  xLrpper)/2.0 

yMid  =  f(xMid) 

iters  =  0 

While  ( (xUpper  -  xLower)/2.0  <  xTol ) 
iters  =  iters  +  1 

if(  yLower  *  yMid  >  0.0)  Then  xLower  =  xMid 

Else  xUpper  =  xMid 

Endofif 

xMid  =  (xEower  +  xUpper)/2.0 
yMid  =  f(xMid) 

Endofwhile 

Return  xMid,  yMid,  iters  (*  xMid  =  approximation  to  root  *) 

Fig.  1  Pseudocode  for  bisection  root  finding  aigorithm 


To  validate  that  the  present  thermokinetic  model  agrees  with  that  of  Chen  et  al.,^^’"^^  the 
thermokinetic  model  was  plotted  using  the  parameters  from  Tables  1  and  2  along  with  the 
experimental  data  from  Darling  et  al.'^  in  Fig.  2.  The  saturated  formulation  of  the  thermokinetic 
model  was  used,  although  the  unsaturated  model  yields  the  same  curve  with  slightly  different 
values  of  AHseg  (+3.5%)  and  Ql  (—3.0%).  The  thermokinetic  model  agrees  well  with  the 
experimental  data  and  the  model  implementation  agrees  with  that  of  Chen  et  al.^^’"^^  Moreover, 
this  model  for  nanocrystalline  grain  growth  assumes  both  kinetic  and  thermal  contributions,  in 
agreement  with  experimental  findings  of  Darling  et  al.^^  In  the  remainder  of  this  report,  we 
examine  the  sensitivity  of  the  thermokinetic  model  to  input  parameters  and  discuss  extending  it 
towards  predicting  nanocrystalline  grain  growth  at  a  variety  of  temperatures  and  solute 
concentrations. 

Optimization  in  the  present  analysis  uses  an  unconstrained  nonlinear  optimization  algorithm. 
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Fig.  2  Evolution  of  the  average  grain  size  as  a  function  of  annealing 
time  for  nanocrystalline  Fe-4%  Zr  alloy  annealed  at  800  °C 
(1,073  K) 


i.e.,  the  Nelder-Mead  simplex  algorithm."*^  This  algorithm  uses  a  simplex  of  n  +  1  points  for 
n-dimensional  veetors  x.  The  algorithm  first  makes  a  simplex  around  the  initial  guess  xq  by 
adding  5%  of  eaeh  eomponent  xo{i)  to  xq,  and  using  these  n  veetors  as  elements  of  the  simplex  in 
addition  to  Xq.  Then,  the  algorithm  modifies  the  simplex  repeatedly  until  it  meets  a  stopping 
eriterion  and  obtains  a  solution  Xgoi-  Some  of  the  basie  MATLAB  eodes  for  the  thermokinetie 
model  and  its  optimization  with  the  temperature  dependenee  are  attaehed  in  Appendiees  A-C. 
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3.  Experimental  Fe-Zr  Results 


3.1  X-ray  Diffraction  (XRD)  Measurements  of  Grain  Growth  in  Nanocrystalline  Fe 

This  section  discusses  a  newly-acquired  XRD  dataset  for  Fe-Zr  alloys  with  various 
concentrations  of  Zr  (0.5%  Zr,  1%  Zr,  4%  Zr,  10%  Zr)  and  at  various  annealing  temperatures 
(700,  800,  900,  and  1,000  °C).  Currently,  the  thermokinetic  model  for  Fe-Zr  has  only  been 
applied  to  a  nanocrystalline  Fe-4%  Zr  alloy  at  800  °C.  To  illustrate  the  difficulty  of 
high-temperature  processing  of  nanocrystalline  materials  in  pure  a-Fe  and  to  establish  a  baseline 
for  grain  growth,  a  nanocrystalline-Fe  sample  with  a  nominal  grain  size  of  5  nm  was  annealed  at 
temperatures  of  400,  500,  600,  and  700  °C  (Fig.  3).  The  lines  are  included  as  a  guide  to  the  eye 
for  the  XRD  data  at  different  annealing  temperatures.  At  temperatures  of  greater  than  400  °C, 
there  is  appreciable  grain  growth  in  the  nanocrystalline  Fe  material;  at  temperatures  of  600  and 
700  °C,  there  is  more  than  an  order  of  magnitude  of  difference  in  grain  size  after  annealing  for  1  h 
(3600  s).  These  results  are  in  good  agreement  with  prior  work  on  nanocrystalline  Fe,  which  has 
exhibited  grain  growth  at  temperatures  around  500  °C.  In  the  following  section,  we  discuss  the 
process  used  for  the  nanocrystalline  structures,  the  microstructure  and  operating  mechanisms  for 
grain  growth,  and  the  data  obtained  for  various  concentrations  of  Zr  in  the  Fe-Zr  alloys. 
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Fig.  3  Evolution  of  the  average  grain  size  as  a  function  of  annealing 
time  at  multiple  temperatures  (400-700  °C)  for  unalloyed 
nanocrystalline  Fe 


3.2  Grain  Growth  in  Nanocrystalline  Fe 

The  processing  and  microstructure  of  these  samples  helps  explain  the  observed  X-ray  behavior. 
During  high-energy  mechanical  alloying,  nanoscale  microstructures  are  generated  by  the 
accumulation  of  severe,  localized  plastic  deformation.  After  reaching  a  certain  level  of  strain, 
dislocation  cell  structures  and  small  angle  GBs  approximately  25  nm  in  size  form.  Continued 
mechanical  attrition  forms  smaller  grains  within  the  entire  volume  of  the  material.'*^’'^^  This  initial 
grain  size  distribution  changes  upon  grain  growth  at  high  temperatures;  the  magnitude  and 
direction  of  the  grain  size  distribution  is  dependent  on  the  evolved  representative  population.  The 
change  in  population  can  occur  in  2  primary  ways:  normal  and  abnormal  grain  growth.^*’ 

An  ideal  grain  structure  exhibits  normal  microstructure  growth  that  occurs  continuously  and 
uniformly;  this  structure  is  also  held  in  this  state  by  a  radius -dependent  counterbalancing  effect. 
However,  grains  with  diameters  much  larger  or  smaller  than  the  normal  size  distribution  grow  at 
rates  that  eventually  cause  them  to  rejoin  the  collective  population. In  this  basic  model  of  grain 
growth,  the  GB  area  reduction  and  the  corresponding  energy  reduction  drives  the  growth  rate.  In 
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nanocrystalline  materials,  size-dependent  microstructure  features,  such  as  triple  lines  and 
quadruple  points,  can  kinetically  hinder  self-similar  grain  growth  as  well.^^  Abnormal  grain 
growth,  however,  proceeds  in  a  much  different  manner.  Here,  the  controlling  mechanisms  allow 
select  grains  some  growth  advantage  over  neighboring  grains;  consequently,  accelerated  growth 
of  these  select  grains  occurs  at  the  cost  of  consuming  the  much  smaller  surrounding  grains. 

Abnormal  grain  growth  plays  a  significant  role  in  controlling  the  stability  and  growth  kinetics  of 
nanocrystalline  microstructures.  In  this  work,  the  grain  growth  observed  in  nanocrystalline  Fe 
(Fig.  3)  can  be  controlled  through  adding  a  solute  element  (e.g.,  Zr)  that  segregates  to  the  GBs. 
Hence,  in  the  Fe-Zr  system  reported  here,  grain  growth  is  controlled  by  solute-influenced 
mechanisms,  such  as  second  phase  particle  pinning  and/or  inhomogeneous  solute  drag  effects. 
Inhomogeneous  solute  drag  effects  can  arise  from  a  high  dependency  of  GB  composition  on  its 
inherent  structure  and  orientation.^^  In  a  sample  composed  of  a  large  distribution  of  GB 
configurations,  local  regions  of  anisotropy  in  GB  composition  exist.  GBs  represent  a  finite 
volume  (e.g.,  based  on  GB  area  and  width)  within  the  material  for  segregating  elements  to 
occupy.  At  small  grains  sizes  (more  volume)  and  low  solute  concentrations  (less  segregant),  the 
distribution  of  the  segregating  species  to  the  boundaries  is  expected  to  be  less  uniform,  which 
leads  to  more  abnormal  grain  growth. 

At  elevated  temperatures,  nanocrystalline  Fe  and  Fe-Zr  alloys  undergo  bimodal  grain  growth 
where  both  larger,  abnormal  grains  and  nanoscale  grains  coexist.  While  transmission  electron 
microscopy  can  supply  detailed  information  about  grain  size  over  a  relatively  small  area,  the 
X-ray  diffraction  technique  used  in  the  present  study  can  supply  in  situ  information  about  the 
volume- averaged  grain  size  and  grain  growth  over  much  larger  areas,  encompassing  both  of  the 
aforementioned  modes  of  grain  growth.  X-ray  line  broadening  analysis  takes  advantage  of  the  fact 
that,  below  approximately  100  nm,  small  crystallites  cause  broadening  or  an  angular  divergence 
of  the  diffracted  beam  at  angles  near  the  specific  Bragg  angle.^"^  Therefore,  in  a  given  sample 
volume,  regions  containing  crystallites  much  larger  100  nm  do  not  contribute  to  the  overall 
broadening  observed.  Thus,  line  broadening  analysis  for  grain  size  estimates  is  sensitive  only  to 
the  volume  fraction  of  the  microstructure  containing  nanocrystalline  grains.  Although,  it  should 
be  noted  that  the  accuracy  of  the  X-ray  diffraction  technique  may  decrease  for  average  grain  sizes 
larger  than  50  nm  as  abnormal  grain  growth  and  second  phase  formation  may  occur. 

Further  supporting  the  XRD  measurements  in  Fig.  3,  recent  results  have  shown  that 
nanocrystalline  Fe,  prepared  by  mechanical  milling  and  annealed  isothermally  in  a  tube  furnace 
under  forming  gas,  undergoes  extensive  growth  abnormally  at  600  The  resultant 
microstructure  was  found  to  consist  of  regions  containing  50%  ultrafine  (100-300  nm)  grains  and 
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50%  coarse,  micron  grains  (approximately  6  pm).  By  700  °C,  nanocrystalline  Fe  was  shown  to  be 
entirely  transformed  to  a  microstructure  containing  micron-sized  grains. Figure  3  reveals  that  in 
situ  heating  during  X-ray  analysis  has  similar  results  to  that  reported  earlier. 

The  microstructure  evolution  for  pure,  nanocrystalline  Fe,  after  annealing  at  600  and  700  °C 
during  high-temperature  X-ray  experiments,  clearly  displays  both  normal  and  abnormal  grain 
growth  (Fig.  4).  These  images  are  etched  to  reveal  regions  of  abnormal  grain  growth,  which  are 
the  white  or  gray  regions  with  smooth  contrast  (indicated  by  white  arrows),  while  the  nanometer 
or  ultrafine  grained  regions  are  the  gray  regions  with  a  rough  or  granulated  contrast  (indicated  by 
black  arrows). Figure  4a,  taken  after  annealing  pure  nanocrystalline  Fe  at  600  °C,  shows  that 
the  microstructure  has  undergone  extensive  grain  growth.  At  600  °C,  abnormally  grown  grains 
(average  size  of  15  pm)  represent  a  significant  fraction  of  the  microstructure  while  the  remaining 
fraction  is  nanocrystalline  or  ultrafine  grained.  Annealing  at  700  °C  (Fig.  4b)  caused  the  entire 
microstructure  to  undergo  abnormal  grain  growth,  resulting  in  an  average  grain  size  much  larger 
than  at  600  °C. 


Fig.  4  Evolution  of  the  microstructure  as  a  function  of  annealing  time  at  multiple 

temperatures  (600,  700,  and  1,000  °C)  for  nanocrystalline  Fe  and  Fe-0.5%  Zr  alloy 
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3.3  XRD  Measurements  of  Grain  Growth  in  Nanocrystalline  Fe-Zr 


The  XRD  datasets  for  nanoerystalline  Fe-Zr  alloys  at  higher  temperatures  than  for 
nanoerystalline  Fe  were  also  obtained.  The  volume- averaged  grain  size  data  as  a  funetion  of 
annealing  time  are  shown  in  Fig.  5  for  the  various  Fe-Zr  alloys  and  temperatures.  It  is  important 
to  point  out  a  few  differenees  between  Fig.  3  and  Fig.  5.  First,  the  seale  bar  for  the  grain  size  is 
exaetly  1  order  of  magnitude  lower  in  Fig.  5  (i.e.,  the  maximum  grain  size  for  the  y-axis  was 
redueed  from  600  nm  to  60  nm).  Seeond,  the  annealing  temperatures  are  mueh  higher  in  Fig.  5 
(i.e.,  the  range  of  temperatures  was  inereased  from  400-700  °C  to  700-1,000  °C).  Interestingly, 
the  smallest  (0.5%  Zr)  and  largest  (10%  Zr)  eoneentrations  of  Zr  have  the  smallest  stabilized 
grain  sizes  for  700-900  °C  range.  However,  in  both  of  these  Fe-Zr  systems,  there  is  very  different 
behavior  at  1,000  °C,  whieh  may  be  attributed  to  seeond  phase  formation  (10%  Zr)  or  suffleient 
kinetie  energy  to  overeome  the  influenee  of  Zr  on  stabilizing  the  GBs  (0.5%  Zr).  The  present 
dataset  is  used  for  fitting  the  thermokinetie  model. 
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Grain  Size  (nm)  Grain  Size  (nm) 


(a)  Fe-0.5 %  Zr  (b)  Fe-1  %  Zr 


(c)  Fe-4%  Zr  (d)  Fe-10%  Zr 

Fig.  5  Evolution  of  the  average  grain  size  as  a  function  of  annealing  time  at  multiple  temperatures 

(700-1,000  °C)  for  nanocrystalline  Fe-Zr  alloy  with  different  Zr  concentrations:  a)  0.5%  Zr,  b) 
1%  Zr,  c)  4%  Zr,  d)  10%  Zr 
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3.4  Grain  Growth  in  Nanocrystalline  Fe-Zr 


The  addition  of  Zr  as  a  solute  in  Fe  has  been  shown  to  stabilize  the  grain  size  to  high  homologous 
temperatures. The  stability  was  observed  to  increase  with  Zr  content,  up  to  10  at%  Zr  where 
the  precipitation  of  equilibrium  intermetallics  was  exceedingly  favorable.'^  The  increase  in 
stability  with  solute  content  followed  by  precipitation  of  intermetallics  and  loss  of  stability  have 
been  used  to  delineate  a  system  and  its  extent  of  being  thermodynamically  stable  by  many 
authors.^*  It  is  thus  expected  that  the  concentration  of  stabilizing  solute  plays  a  critical  role  in  the 
overall  evolution  of  such  microstructures  at  elevated  temperatures.  Darling  et  al.  have  reported 
that  an  increase  in  Zr  content  leads  to  a  smaller  retained  grain  size  within  the  nanoscale  regions  of 
the  microstructure  and  to  some  extent  controls  the  level  of  abnormal  grain  growth.  Additionally, 
the  concentration  of  the  Zr  solute  has  been  observed  to  influence  the  grain  morphology.  Low  Zr 
contents  (1/3  at%)  led  to  acicular-type  grains  with  aspect  ratios  greater  than  5  in  some  cases, 
while  higher  Zr  concentrations  (1-4  at%)  led  to  more  equiaxed  structures.  The  amount  of  solute 
also  directly  impacts  the  as-milled  grain  size.^^’^^ 

In  addition  to  compositional  effects,  it  has  also  been  proposed  that  the  body-centered  cubic  (bee) 
to  face-centered  cubic  (fee)  phase  transition  may  be  another  point  of  system  instability. 

Figure  4  is  in  line  with  these  results.  Figure  4c  shows  the  etched  micrographs  of  Fe-0.5  at%  Zr 
after  annealing  at  700  °C  and  1,000  °C,  in  which  there  are  only  a  few  small  regions  that  have 
evidence  of  abnormal  grain  growth.  The  fraction  of  the  microstructure  abnormally  grown  after 
annealing  at  700  °C,  as  estimated  from  lower  magnification  micrographs,  was  found  to  be  less 
than  3%.  This  is  a  stark  improvement  over  pure  Fe.  However,  after  annealing  at  1,000  °C,  the 
Fe-0.5  at%  Zr  alloy  (Fig.  4d)  has  undergone  extensive  abnormal  grain  growth.  Samples 
containing  greater  than  0.5  at%  Zr  show  limited  abnormal  growth  (less  than  5%  of  the  area),  even 
at  temperatures  as  high  as  1,000  °C.  It  is  expected  that  the  10  at%  Zr  sample  is  heavily 
precipitated.  In  this  case,  kinetic  pinning  by  a  secondary  phase  dominates  the  microstructural 
stability;  such  a  mechanism  is  not  captured  by  the  current  model. 

Since  lower  Zr  solute  contents  (0.5  at%)  result  in  bimodal  microstructures  having  large  areas  of 
abnormal  growth  and  regions  where  the  morphology  of  the  nanoscale  grains  may  be  largely 
distorted,  the  reported  model  is  not  expected  to  fit  this  data  well.  Additionally,  the  10  at%  Zr  data 
are  most  likely  heavily  precipitated,  giving  rise  to  Zener  pinning;  such  growth  kinetics  are  also 
not  captured  within  the  current  model.  Hence,  the  XRD  data  from  1  at%  and  4  at%  Zr — where  the 
majority  of  the  grain  growth  is  expected  to  be  normal,  that  is,  controlled  by  both  thermodynamics 
and  kinetics  (solute  drag) — are  ideal  experimental  datasets  for  the  present  model.  It  is  expected 
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that  the  1  at%  and  4  at%  Zr  samples  contain  sufficient  solute  to  evenly  coat  the  GB,  without  being 
high  enough  to  kinetically  promote  precipitation  of  secondary  phases. 


4.  Model  Results 


4.1  Parameter  Sensitivity  Analysis 

The  sensitivity  analysis  of  a  model  is  a  local  measure  of  the  effect  of  a  given  input  (or  parameter) 
on  a  given  output  (in  this  case,  the  evolution  of  grain  size  with  time).  Another  definition  of 
sensitivity  analysis  is,  “The  study  of  how  the  uncertainty  in  the  output  of  a  model  can  be 
apportioned  to  difference  sources  of  uncertainty  in  the  model  input”,  i.e.,  it  is  important  to  know 
which  parameters  affect  the  output  and  the  magnitude  of  this  interaction.  Prior  to  fitting  to  the 
XRD  datasets  described  in  Section  3,  some  of  the  material  parameters,  which  were  used  within  the 
thermokinetic  model  for  Fe-4%  Zr  model  at  800  °C,  were  perturbed  to  understand  their  influence 
on  the  grain  size  and  the  shape  of  the  grain  size-annealing  time  curve.  There  are  several  methods 
for  examining  parameter  sensitivity  since  the  sensitivity  is  typically  defined  as  Sj  =  dY/dXj, 
where  Y  is  the  output  and  Xj  is  the  jth  parameter.  From  a  local  perspective,  parameter  Xj  can  be 
locally  adjusted  around  some  default  value  while  all  other  parameters  are  held  fixed  to  quantify 
the  change  in  output  Y.  From  a  global  perspective,  though,  it  may  be  important  to  know  the  effect 
of  parameter  Xj  on  the  output  Y  while  all  other  parameters  Xi^j  are  also  varying  (e.g.,  a  Monte 
Carlo  approach  or  the  Morris  one-factor-at-a-time  approach).  It  should  be  noted  that  instead  of  a 
scalar  value  for  the  output  Y,  in  this  study  F  is  a  function  of  time;  hence,  the  sensitivity  of  grain 
size  to  a  particular  parameter  may  change  with  time.  The  local  and  global  sensitivity  analyses  are 
briefly  explored  in  Sections  4.1.1  and  4.1.2.  Additionally,  the  change  in  this  behavior  as  a 
function  of  different  temperatures  and  concentrations  of  Zr  is  also  investigated  in  Section  4.1.3. 

4.1.1  Local  Sensitivity  Analysis 

The  local  sensitivity  of  the  parameters  was  examined  by  selecting  a  baseline  set  of  parameter 
values  (the  Fe-4%  Zr  model  at  800  °C)  and  iteratively  varying  each  by  different  percentages  (e.g., 
-1%,  -1-1%,  -1-2%,  -1-3%)  to  qualitatively  observe  their  effect  on  the  grain  size  evolution  curves.  The 
subsequent  global  sensitivity  analysis  in  Section  4.1.2  will  use  a  Monte  Carlo  analysis  to  quantify 
the  sensitivity  of  grain  size  to  each  parameter  (in  combination  with  the  other  parameters  varying) 
as  a  function  of  time. 

The  change  in  the  grain  growth  behavior  as  a  function  of  deviations  from  the  material  parameters 
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utilized  in  the  present  Fe-4%  Zr  model  at  800  °C  (i.e.,  70,  AHgeg,  Tm,  and  Ql)  are  shown  in 
Fig.  6.  In  each  subplot,  the  blue  line  represents  the  baseline  and  the  green  and  red  lines  indicate 
deviations  from  these  baseline  values  (Table  2)  in  increments  of  +1%  and  —1%,  respectively. 
First,  increasing  70  results  in  an  increase  in  the  equilibrium  grain  size  D*  and  a  corresponding 
increase  in  grain  growth  behavior.  However,  for  all  other  parameters,  an  increase  corresponds  to  a 
decrease  in  both  D*  and  the  grain  growth  behavior.  The  time  scale  associated  with  reaching  the 
stabilized  grain  size  correspondingly  decreases/increases  with  D*.  Second,  the  rate  of  deviation 
from  the  baseline  grain  growth  plot  (in  terms  of  deviation  at  large  times)  for  a  1%  deviation  in 
material  parameters  (i.e.,  the  sensitivity)  follows  this  order:  AHgeg  (largest),  70,  Fm,  and  Ql 
(smallest).  However,  while  the  activation  energy  Ql  has  the  smallest  effect  at  larger  times,  it  has  a 
dramatic  effect  on  the  shape  of  the  grain  growth  curve  at  shorter  times  (i.e.,  there  is  a  greater 
sensitivity  to  Ql  at  shorter  times).  In  fact,  for  lower  activation  energies,  the  thermokinetic  model 
starts  to  approach  the  behavior  of  the  thermokinetic  model  without  the  GB  energy  effect,  i.e.,  the 
Gong  et  al.model.^^  Notice  that,  in  some  cases,  only  a  few  red  or  green  curves  are  generated.  In 
each  case,  10  increments  of  ±1%  were  attempted.  However,  there  is  a  limit  that  these  variables 
can  be  increased/decreased  that  corresponds  to  a  regime  where  the  stabilized  grain  size 
D*  (=  —72/71)  becomes  either  lower  than  the  initial  grain  size  (i.e.,  72  <  —Dq'Ji,  no  grain 
growth)  or  negative  (i.e.,  both  71  and  72  are  positive,  no  stabilized  grain  size). 
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Grain  Size  (nm)  Grain  Size  (nm) 


(a)  7o  (b)  I^Hseg 


(c)  Tfeo 


(d)  Ql 


Fig.  6  Evolution  of  the  average  grain  size  as  a  function  of  annealing  time  for  nanocrystalline  Fe^%  Zr 
alloy  annealed  at  800  °C  (1,073  K).  The  different  plots  (a-d)  show  the  thermokinetic  model  with 
the  parameters  listed  in  Tables  1  and  2  and  the  new  curves  generated  by  changing  the  value  of  one 
of  the  parameters  (in  increments  of  ±1%):  a)  70,  b)  AHgeg,  c)  Tbo,  and  d)  Ql. 
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The  rate  of  grain  growth  ean  also  be  ealeulated  from  the  preeeding  analysis  as  a  funetion  of 
annealing  time,  as  shown  in  Fig.  7.  A  eentral  differenee  numerieal  method  was  used  to  ealeulate 
the  derivative  of  the  grain  growth  eurves  (Fig.  6).  Similar  trends  are  observed.  Interestingly,  the 
aetivation  energy  again  shows  a  large  influenee  on  the  rate  of  grain  growth  in  the  beginning  stages 
of  annealing.  In  faet,  while  AHs^g,  70,  and  F^o  show  a  sealing  behavior  in  terms  of  the  rate  of 
grain  growth,  the  Ql  plot  shows  very  different  behavior.  For  inereasing  Ql,  the  rate  of  grain 
growth  is  deereased  as  with  the  other  material  parameters.  However,  for  deereasing  Ql,  there  are 
2  regimes:  a  regime  of  aeeelerated  grain  growth  rates  at  short  times,  and  a  regime  where  grain 
growth  is  lower  than  the  baseline.  In  this  manner,  the  aetivation  energy  Q l  controls  the  kinetics 
for  obtaining  the  stabilized  grain  size  without  actually  modifying  D*  (i.e.,  modifies  the  shape  of 
the  grain  growth  curve,  as  can  be  seen  in  Fig.  6). 
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Rate  of  Grain  Growth  (nm/s)  Rate  of  Grain  Growth  (nm/s) 


(a)  7o  (b)  I^Hseg 


(c)  Tfeo  (d)  Ql 

Fig.  7  Evolution  of  the  rate  of  grain  growth  as  a  function  of  annealing  time  for  nanocrystalline  Fe-4%  Zr 
alloy  annealed  at  800  °C  (1,073  K).  In  a  similar  manner  to  Fig.  6,  the  different  plots  (a-d)  show  the 
thermokinetic  model  with  the  parameters  listed  in  Tables  1  and  2  and  the  new  curves  generated  by 
changing  the  value  of  one  of  the  parameters  (in  increments  of  ±1%):  a)  70,  b)  AHseg,  c)  Fbo,  and 
d)  Ql. 
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Additionally,  the  influence  of  additional  parameters  in  the  thermokinetic  model  can  be  explored. 
For  instance,  the  effect  of  the  mobility  term,  Mq,  on  the  grain  growth  curves  is  shown  in  Fig.  8.  In 
these  plots,  10  increments  of  ±10%  were  used  to  explore  the  effect  of  the  mobility  term  on  grain 
growth.  Hence,  the  upper  bound  for  the  green  curves  represents  doubling  the  mobility  term 
(±100%)  and  the  lower  bound  of  the  red  curves  represents  a  mobility  of  0  (—100%),  i.e.,  no  grain 
growth.  Interestingly,  the  mobility  behaves  very  similar  to  the  activation  energy  Ql;  high 
mobilities  favor  fast  grain  growth  without  changing  the  stabilized  grain  size.  However,  while  the 
mobility  results  in  changes  similar  to  the  activation  energy,  much  larger  increases  in  mobility  are 
required  to  make  a  substantial  change  in  the  grain  growth  behavior  (i.e.,  the  thermokinetic  model 
is  relatively  insensitive  to  Mq). 


(a)  (b) 

Fig.  8  Evolution  of  the  average  grain  size  as  a  function  of  annealing  time  for  nanocrystalline  Fe^%  Zr 
alloy  annealed  at  800  °C  (1,073  K).  The  different  plots  show  the  thermokinetic  model  with  the 
parameters  listed  in  Tables  1  and  2  and  the  new  curves  generated  by  changing  the  mobility  term 
Mq  (in  10  increments  of  ±10%). 


The  preceding  analysis  shows  that  a  fit  of  the  material  parameters  to  experimental  grain  growth 
data  does  not  result  in  a  unique  set  of  material  parameters  (e.g.,  70,  T^q,  AHs^g,  Ql).  For  instance, 
using  the  parameters  outlined  for  the  Fe-4%  Zr  at  800  °C,  the  GB  energy  term  70  was  modified 
by  multiplying  it  by  0.97  to  1.03  in  increments  of  0.01.  The  initial  curve  was  used  as  a  baseline 
curve  and  a  curve  with  the  new  value  of  70  was  fit  to  the  original  curve  using  only  the  enthalpy 
change  AHgeg  and  activation  energy  Ql  (Table  3).  The  new  curves  were  virtually 


22 


indistinguishable  from  the  baseline  grain  growth  curve.  Interestingly,  the  activation  energy  Ql 
was  unchanged  and  only  the  enthalpy  change  /S.Hseg  was  modified.  This  result  agrees  with  Fig.  6; 
i.e.,  the  activation  energy  modifies  the  shape  of  the  curve  and  70,  AHseg,  and  Fbo  all  modify  the 
magnitude  of  the  grain  size  for  the  thermokinetic  model.  Moreover,  there  is  a  linear  dependence 
between  70  and  AHs^g  for  non-unique  solutions,  e.g.,  70  =  cAHgeg  with  c  = 

8.175  X  10“®  molm“^.  The  saturated  excess  Tm  was  held  constant  in  the  present  analysis,  but  Tm 
is  likely  to  have  a  similar  effect  to  AHseg  based  on  Fig.  6. 


Table  3  Non-unique  material  parameter  set  for  Fe^%  Zr  at  800  °C 


Parameter 

Symbol 

0.9770 

0.9870 

0.9970 

70 

I.OI70 

1.0270 

1.0370 

GB  energy  (J  m“^) 

7o 

0.708 

0.715 

0.723 

0.73 

0.737 

0.745 

0.752 

Saturated  excess  (10“^  molm“^) 

Ffco 

1.25 

1.25 

1.25 

1.25 

1.25 

1.25 

1.25 

Enthalpy  change  ( kJ  mol~^) 

AHseg 

87.5 

88.1 

88.7 

89.3 

89.9 

90.5 

91.1 

Activation  energy  (kJ  mol”^) 

Ql 

345.1 

345.1 

345.1 

345.1 

345.1 

345.1 

345.1 

4.1.2  Global  Sensitivity  Analysis 

The  global  sensitivity  analysis  of  the  parameters  can  be  explored  using  a  number  of  different 
methods.  Herein,  we  have  used  a  Monte  Carlo  approach  that  has  been  coupled  with  scatter  plots, 
regression  modeling,  and  correlation  coefficients  to  show  the  sensitivity  of  the  grain  size 
evolution  to  the  model  parameters.  For  the  present  analysis,  6  parameters  were  varied  within  the 
Monte  Carlo  simulations:  AHseg  (heat  of  segregation),  Ql  (activation  energy),  70  (initial  grain 
boundary  energy),  F^o  (saturated  solute  excess),  Mq  (grain  boundary  mobility),  and  Dq  (initial 
grain  size).  The  last  2  parameters  were  included  to  show  the  relative  sensitivity  of  the 
thermodynamic  and  kinetic  parameters  (AHseg,  Ql,  lo,  Tfeo)  to  the  material-based  parameters 
(Mq,  Dq).  The  computationally-inexpensive  thermokinetic  grain  growth  model  allowed  for  a  large 
sample  size  (N  =10,000  simulations)  to  be  used  with  parameter  values  chosen  randomly  within 
the  bounds  of  -0.5%  and  -i-0.5%  of  the  default  values.  Varying  within  a  set  range  of  percentages 
essentially  normalizes  each  parameter  so  that  their  sensitivities  can  be  prepared.  A  summary  of 
the  grain  growth  curves  is  shown  in  Fig.  9,  with  the  2%,  10%,  and  50%  percentile  (median) 
curves  shown  as  a  guide  to  the  shape  of  the  distribution.  These  curves  can  then  be  analyzed  for 
their  sensitivity  with  respect  to  the  various  Monte  Carlo  parameters. 
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Fig.  9  Grain  size  evolution  as  a  function  of  annealing  time  for  the  Monte  Carlo 
analysis 


Scatter  plots  are  occasionally  used  to  visually  examine  the  correlation  between  a  model  parameter 
and  the  output.  For  instance,  Fig.  10  shows  the  scatter  in  grain  size  at  1  h  as  a  function  of  the  4 
thermodynamic  and  kinetic  parameters.  The  red  markers  represent  the  10,000  data  points  and  the 
line  represents  the  best  fit  linear  regression  model  for  each  of  the  different  parameters  (without 
including  the  other  parameters).  The  spread  of  the  data  points  indicates  the  variation  as  a  function 
of  the  other  parameters  —  the  spread/variation  is  wider  for  some  of  the  parameters  because  they 
do  not  account  for  the  variation  in  the  larger  effects  (higher  sensitivity  parameters).  Another 
method  for  examining  the  correlation  between  parameters  and  responses  is  to  perform  a  rank 
transformation;  in  this  respect,  any  nonlinear  differences  because  of  the  magnitude  of  the 
parameter  or  response  changes  is  removed  such  that  the  monotonic  nature  of  the  relationships  are 
examined.  The  rank  transformed  scatter  plot  is  shown  in  Fig.  1 1.  In  Figs.  10  and  11,  Hgeg  has  the 
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largest  effect  (based  on  slope  of  the  linear  regression  line),  followed  by  70  and  Fbo,  with  very  little 
effect  from  Ql.  Again,  the  largest  effect  of  Ql  is  at  shorter  times  than  the  1  h  time  chosen  here, 
so  this  effect  isn’t  surprising.  The  rank  transformation  shows  much  wider  scatter  than  the 
untransformed  data. 


(c)  (d) 

Fig.  10  Scatter  plot  of  the  grain  size  as  a  function  of  various  thermodynamic  and  kinetic  parameters:  a) 
Hseg,  b)  Ql,  c)  7o,  and  d)  r^o-  The  scatter  of  the  data  points  is  due  to  the  variation  in  the  other 
parameters  at  1  h.  The  line  represents  a  best  fit  linear  regression  model. 
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(a)  (b) 


(c)  (d) 

Fig.  1 1  Rank  transformation  scatter  plot  of  the  grain  size  as  a  function  of  thermodynamic  and  kinetic 

parameters:  a)  Hseg,  b)  Ql,  c)  70,  and  d)  r^o-  The  scatter  of  the  data  points  is  due  to  the  variation 
in  the  other  parameters  at  1  h.  The  line  represents  a  best  fit  linear  regression  model  to  the  rank 
transformation  data. 
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In  Fig.  12,  the  correlation  coefficient  (R)  evolution  for  each  of  the  parameters  is  plotted  as  a 
function  of  time.  This  method  gives  both  the  direction  of  the  effect  and  the  magnitude  of  the 
effect  of  each  parameter.  In  addition  to  R,  the  p- value  was  also  calculated.  The  parameters  Mq 
and  Do  were  not  statistically  significant  (based  on  a  p- value  of  0.05),  the  activation  energy  Ql 
was  only  statistically  significant  for  short  times  (t  <f»2  h),  and  the  remaining  parameters  were 
statistically  significant  for  all  t. 
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Fig.  12  Evolution  of  the  correlation  coefficients  as  a  function  of  annealing  time  for 
the  Monte  Carlo  analysis 


The  global  sensitivity  of  the  6  parameters  can  be  assessed  by  fitting  the  10,000  grain  size 
evolution  curves  to  the  parameters.  In  this  respect,  at  any  time,  the  output  (grain  size)  Y  can  be 
modeled  as  a  function  of  the  parameters  x;  i.e.,  Y  =  f  (x.)  =  f  {xi,  X2,  ■  ■  ■ ,  Xj)  where  j  is  the 
number  of  parameters.  In  this  respect,  the  grain  size  can  be  modeled  by  the  linear  regression 
model: 


Y  =  Co  +  ^a 


i=l 


^  _  1 


(17) 


whereby  Cq  is  the  average  grain  size,  Ci  is  the  coefficient  (slope)  for  parameter  i,  and  the  final 
term  in  the  summation  normalizes  the  parameters  of  x,  by  their  default  value  (and  subtracts  1)  to 
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give  a  percentage  increase  or  decrease  (i.e.,  normalizing  each  parameter  to  be  able  to  compare 
their  effect  on  grain  size  Ci).  The  coefficient  C*  is  similar  to  the  sensitivity  since  it  is  the 
derivative  of  the  grain  size  with  respect  to  the  percentage  increase/decrease.  Here,  we  differentiate 
the  two  by  taking  the  absolute  value  of  C*  for  a  measure  of  the  sensitivity,  i.e.,  abs  (Q).  In  this 
respect,  Ci  shows  the  direction  of  the  parameter  effect  while  abs  {Ci)  shows  the  degree  of 
sensitivity  for  comparison  of  the  different  parameters. 

The  evolution  of  the  regression  model  coefficients  is  plotted  as  a  function  of  time  to  show  the 
sensitivity  of  the  grain  size  evolution  to  these  parameters.  Figure  13  shows  the  evolution  of  the 
regression  model  coefficients  Ci  (and  the  parameter  sensitivity)  as  a  function  of  time.  First, 

Fig.  13a  shows  the  evolution  of  C*  for  each  of  the  parameters  as  a  function  of  time.  The  grain 
boundary  energy  70  is  the  only  parameter  that  causes  the  entire  grain  size  evolution  curve  to 
increase  with  an  increase  to  the  parameter  (i.e.,  there  is  a  direct  correlation  with  70).  Additionally, 
the  mobility  Mq  and  initial  grain  size  Dq  have  almost  no  effect  and  the  activation  energy  Q l  only 
has  an  effect  at  short  times.  In  Fig.  I3b,  abs  (C*)  is  plotted  for  the  6  parameters.  At  very  short 
times,  the  activation  energy  Q l  can  have  as  much  of  an  effect  as  some  of  the  other 
thermodynamic  and  kinetic  parameters,  but  this  effect  decreases  as  time  increases.  Moreover,  the 
sensitivity  of  the  heat  of  segregation  Hgeg  is  the  largest,  followed  by  the  grain  boundary  energy  (Tq 
and  the  saturated  solute  excess  Ffeo-  However,  the  evolution  of  the  sensitivity  with  time  is  similar 
for  Hseg,  fTo,  and  Ffeo,  as  is  apparent  when  these  parameters  are  normalized  by  their  maximum 
sensitivity  value  (Fig.  13c)  and  when  the  evolution  of  the  normalized  model  coefficient  Ci  for 
Hseg  is  compared  against  that  of  the  other  parameters  (Fig.  13d).  Hence,  changes  in  the 
parameters  Hseg,  ctq,  and  /Fm  result  in  a  nearly  identical  effect  on  the  grain  size  evolution. 

Clearly,  the  global  sensitivity  analysis  of  Fig.  13  shows  that  only  2  parameters  have  independent 
effects  on  the  thermokinetic  model  (Ql  and  Hseg/cro/Tbo)- 


28 


Normalized  Absolute  Model  Coefficient,  C.  Coefficient,  C. 


(a) 


(b) 


(c) 


(d) 


Fig.  13  Grain  size  sensitivity  as  a  function  of  annealing  time  for  the  Monte  Carlo  analysis  of  the  different 
model  parameters 


29 


There  are  a  few  caveats  about  the  present  analysis  that  should  be  mentioned.  First,  the  6 
parameters  explored  were  each  varied  by  the  same  percent  difference  to  explore  their  sensitivity 
with  respect  to  this  difference.  In  many  cases,  if  the  range  and  distribution  of  a  particular 
parameter  is  known  and  can  be  measured,  this  uncertainty  can  instead  be  used  for  the  Monte 
Carlo  analysis.  As  an  example,  as  shown  in  Fig.  8,  certainly  Mq  has  an  effect  if  it  is  doubled,  but 
it  has  no  statistical  significance  if  modified  by  ±0.5%  with  all  other  parameters  varying  by  the 
same  percentage.  Second,  other  sensitivity  analysis  methods  were  also  explored  (e.g.,  the  Morris 
one-at-a-time  method)  and  resulted  in  similar  findings.  Furthermore,  in  varying  many  of  the  other 
material  parameters,  their  effect  was  also  minimal  in  comparison  to  the  thermodynamic  and 
kinetic  parameters  of  the  thermokinetic  model. 

4.1.3  State  Variables:  Temperature  and  Concentration 

The  change  in  the  grain  growth  behavior  as  a  function  of  temperature  and  concentration  is  of 
interest  for  ultimately  tailoring  the  properties  of  the  nanocrystalline  Fe-Zr  system.  The  evolution 
of  the  average  grain  size  as  a  function  of  different  temperatures  and  concentrations  is  shown  in 
Fig.  14.  Again,  the  baseline  model  for  nanocrystalline  Fe-4%  Zr  alloy  annealed  at  800  °C 
(1,073  K)  was  used  with  deviations  up  to  50  °C  and  0.5%  Zr  concentration.  The  blue  lines 
represent  a  lower  temperature/concentration  (750  °C  and  3.5%  Zr)  and  the  red  lines  represent  the 
higher  temperature/concentration  (850  °C  and  4.5%  Zr)  with  intermediate  lines  shown  as  well 
(black  dotted  lines).  Lower  temperatures  and  higher  solute  concentrations  have  a  smaller 
stabilized  grain  size  and  reach  this  stabilized  grain  size  much  quicker  than  higher  temperatures 
and  lower  solute  concentrations,  as  expected.  For  example,  higher  solute  concentrations  at  the 
same  temperature  means  that  the  GBs  can  more  quickly  saturate  with  solute,  thus  providing  an 
impediment  to  further  grain  growth.  While  the  trends  seem  appropriate,  it  should  be  noted  that  the 
magnitude  of  these  changes  may  not  be  correct,  though.  Further  data  would  be  needed  to  validate 
this  behavior  as  a  function  of  temperature  and  concentration. 
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(a)  Temperature,  T  (b)  Coneentration,  xq 

Fig.  14  Evolution  of  the  average  grain  size  as  a  function  of  annealing  time  for  a)  nanocrystalline 
Fe-4%  Zr  alloy  annealed  at  various  temperatures  and  b)  nanocrystalline  Fe-Zr  alloy  with 
different  Zr  concentrations  annealed  at  800  °C  (1,073  K) 


The  thermokinetie  model  ean  use  these  state  variables  to  show  potential  proeessing  paths  that  ean 
be  used  for  tailoring  the  grain  size  (and  properties)  of  nanoerystalline  materials.  We  developed  a 
thermokinetie  stability  map  as  a  funetion  of  temperature  T  and  solute  eoneentration  xq  (Fig.  15). 
The  thermokinetie  stability  map  fits  the  present  model  to  the  experimental  data  obtained  (Fig.  2) 
and  ean  then  ealeulate  the  ehange  in  D*  with  temperature  and  eoneentration.  In  eertain  regions, 
the  present  model  prediets  no  grain  growth  and,  in  other  regions,  the  model  prediets  no  stabilized 
grain  size.  However,  it  should  be  noted  that  the  present  model  does  not  aeeount  for  ehanges  in  the 
material  parameters  as  a  funetion  of  temperature  and  eoneentration.  Present  applieations  of  the 
thermokinetie  model  have  been  to  fit  experimental  datasets  in  both  nanoerystalline  metal 
(e.g.,  Fe-Zr,  Ru-Al)  and  eeramie  systems.  In  truth,  extending  these  models  to  inelude  the  ehange 
in  material  parameters,  sueh  as  GB  energy,  with  temperature  and/or  eoneentration  would  be 
benefieial  for  generating  thermokinetie  stability  maps  for  proeessing  nanoerystalline  materials. 
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Fig.  15  Thermokinetic  stability  map  of  the  evolution  of  the  stabilized 
grain  size  as  a  function  of  temperature  T  and  solute 
concentration  xq  for  the  nanocrystalline  Fe-Zr  system.  The 
white  square  data  point  shows  the  location  of  the  experimental 
data  used  for  fitting  the  thermokinetic  model  (in  Fig.  2)  for  a 
nanocrystalline  Fe-4%  Zr  alloy  annealed  at  800  °C  (1,073  K). 


4.2  Thermokinetic  Model  Extension 

This  subsection  focuses  on  extending  the  present  model  to  handle  a  wider  range  of  grain  growth  in 
nanocrystalline  materials.  First,  the  present  model  is  applied  to  the  Fe-4%  Zr  alloy  at  the  multiple 
temperatures  associated  with  the  XRD  data.  The  optimal  material  parameters  are  listed  in  Table  4 
for  3  different  cases:  1)  all  4  parameters  included  in  optimization,  2)  70  is  held  constant,  and  3) 
both  7o  and  F m  are  held  constant.  These  3  cases  are  separated  into  blocks  within  the  table.  Certain 
trends  emerge  based  on  these  optimal  parameter  sets.  In  the  first  case,  it  can  be  seen  that  the  GB 
energy  70  is  strongly  tied  to  the  enthalpy  change  AHgeg-  Hence,  in  the  second  and  third  cases,  the 
GB  energy  is  held  constant  to  examine  the  influence  of  the  other  parameters.  Furthermore,  the 
saturated  excess  F was  trending  in  opposite  directions  with  increasing  temperature  for  the  first  2 
cases  so  that,  in  the  third  case,  this  parameter  was  held  constant.  Interestingly,  in  all  3  cases,  the 
activation  energy  Ql  tends  to  decrease  with  increasing  temperature.  Furthermore,  when  the  GB 
energy  is  held  constant,  the  enthalpy  change  decreases.  Because  the  range  of  applicability  of  these 
models  is  limited  (Fig.  14),  including  the  influence  of  temperature  on  these  material  parameters 
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extends  the  range  of  applicability  for  this  model.  Figure  16  shows  examples  of  the  present 
thermokinetic  model,  without  any  extensions,  fit  to  the  Fe-Zr  grain  size  data  herein. 


Table  4  Optimal  material  parameter  sets  for  Fe-4%  Zr 


Temperature 

GB  energy 

Saturated  exeess 

Enthalpy  ehange 

Aetivation  energy 

CO 

(Jm-2) 

(10“®molm-2) 

(kJmol-i) 

(kJmol-i) 

T 

70 

Fbo 

AHseg 

Ql 

700 

0.729 

1.243 

89.7 

337.5 

800 

0.733 

1.228 

90.6 

332.2 

900 

0.755 

1.159 

95.7 

324.2 

1,000 

0.743 

1.178 

93.6 

327.0 

700 

0.730“ 

1.257 

89.2 

337.7 

800 

0.730“ 

1.247 

89.6 

332.5 

900 

0.730“ 

1.358 

84.6 

327.1 

1,000 

0.730“ 

1.488 

80.1 

331.5 

700 

0.730“ 

1.250“ 

89.5 

337.5 

800 

0.730“ 

1.250“ 

89.4 

332.5 

900 

0.730“ 

1.250“ 

89.1 

325.6 

1,000 

0.730“ 

1.250“ 

89.1 

328.2 

^  Indieates  that  this  variable  was  held  eonstant  during  optimization. 
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Fig.  16  Evolution  of  the  average  grain  size  as  a  function  of  annealing 
time  at  multiple  temperatures  (700-1,000  °C)  for 
nanocrystalline  Fe-4%  Zr.  The  lines  represent  the  fit  from  the 
thermokinetic  model  for  the  various  temperatures. 


4.2.1  Temperature  Dependence 

The  thermokinetic  model  was  extended  by  including  the  dependence  of  the  GB  energy  70  and 
activation  energy  Ql  on  temperature  (Note  that  the  parameter  70  was  chosen  instead  of  AHgeg, 
because  of  its  similarity  to  Ql  compare  Eqs.  6  and  12).  A  linear  dependence  on  temperature  was 
assumed  for  the  evolution  of  GB  energy  70,  which  seems  reasonable  given  recently  simulated 
trends  in  the  GB  free  energy  with  temperature,^^  i.e., 

7o(^)  =  70,0  +  7o,i^  (18) 

where  70,0  is  the  intercept  and  707  is  the  slope  of  70  with  respect  to  temperature.  An  Arrhenius 
relation  was  chosen  to  represent  the  temperature  dependence  of  the  activation  energy  Ql,  i.e., 

0l(r)  =  0L,„exp(^-^)  (19) 
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respectively,  where  (5l,o  is  the  pre-exponential  factor  and  Q l^i  is  an  exponential  constant, 
respectively.  The  thermokinetic  model  with  temperature-dependent  70  and  Q  l  was  applied  to  the 
Fe-4%  Zr  XRD  data  and  the  resulting  fit  is  shown  in  Fig.  17.  The  resulting  parameters  for  this 
model  are  listed  in  Table  5  and  the  calculated  values  of  the  material  parameters  are  listed  in 
Table  6.  While  the  overall  fit  of  this  model  is  less  accurate  than  individually  fitting  each 
temperature,  the  general  agreement  with  the  experiments  is  good  and  the  model  captures  the  trend 
of  the  experimental  data.  Furthermore,  because  70  and  Q l  now  change  as  a  function  of 
temperature  and  the  fit  uses  the  entire  Fe-4%  Zr  XRD  dataset,  this  allows  the  model  to  be 
extended  over  a  much  larger  temperature  range  than  the  thermokinetic  model  as  applied  to 
experimental  data  from  1  temperature.  Additionally,  because  of  the  non-uniqueness  of  solutions 
as  applied  to  1  temperature,  applying  the  models  with  temperature  dependence  to  the  entire  range 
of  data  can  help  to  obtain  unique  sets  of  optimal  material  parameters.  Last,  this  model  was  also 
applied  to  the  Fe-1%  Zr  case  (Fig.  18)  and  shows  reasonable  agreement  with  the  experimental 
XRD  data. 
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Fig.  17  Evolution  of  the  average  grain  size  as  a  function  of  annealing 
time  at  multiple  temperatures  (700-1,000  °C)  for 
nanocrystalline  Fe-4%  Zr.  The  lines  represent  the  fit  from  the 
thermokinetic  model  with  temperature  dependence  for  70  and 
Ql  (Eqs.  18  and  19). 
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Table  5  Optimal  temperature-dependent  eonstants  for 
7o  and  Ql  in  modified  thermokinetie  model 
for  Fe-4%  Zr 


7o,o 

70,1 

Ql,o 

Ql,i 

(Jm-2) 

(Jm-^K-i) 

(kJmol“^) 

(kJmol-i) 

1.227 

-3.39  X  10-^ 

805 

7.85 

Table  6  Optimal  material  parameter  sets  for  Fe^%  Zr  in  temperature-dependent 
thermokinetie  model  (ealculated  using  parameters  in  Table  5) 


Temperature 

(°C) 

T 

GB  energy 
(Jm-2) 

70  (T) 

Saturated  excess 

(10“®molm-2) 

Tfoo 

Enthalpy  change 
(kJmol-i) 

^Hseg 

Activation  energy 
(kJ  mol“^) 
Ql{T) 

700 

0.897 

1.250^ 

100^ 

305.1 

800 

0.863 

1.250^ 

100^ 

333.9 

900 

0.829 

1.250^ 

100^ 

359.9 

1,000 

0.795 

1.250^ 

100^ 

383.4 

^  Indicates  that  this  variable  was  held  constant  during  optimization. 
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Fig.  18  Evolution  of  the  average  grain  size  as  a  function  of  annealing 
time  at  multiple  temperatures  (700-1,000  °C)  for 
nanocrystalline  Fe-1%  Zr.  The  lines  represent  the  fit  from  the 
thermokinetic  model  with  temperature  dependence  for  70  and 
Ql  (Eqs.  18  and  19). 


Figures  17  and  18  experimentally  show  that,  irrespective  of  Zr  concentration,  the  average  grain 
size  within  both  samples  (1  and  4  at%  Zr)  increases  as  the  annealing  temperature  is  increased. 
However,  the  average  grain  size  for  samples  containing  Fe-4  at%  Zr  are  smaller  than  in  the 
Fe-1  at%  Zr  samples  at  equivalent  temperatures.  When  comparing  the  2  sets  of  data  in 
Figs.  17  and  18,  the  experimental  XRD  datasets  reveal  changes  in  the  kinetics  of  grain  growth. 
The  alloy  containing  4  at%  Zr  shows  similar  growth  rates  for  all  temperatures  tested.  A  steady 
increase  in  the  average  grain  size  occurs  for  times  up  to  1,000  s,  after  which  saturation  occurs. 
The  alloy  containing  1  at%  Zr  behaves  similarly  at  700  °C,  but  growth  above  that  temperature  is 
nearly  instantaneous,  reaching  saturation  within  the  first  30  s. 

4.3  Grain  Size  Distribution  Dependence 

It  is  also  important  to  understand  how  an  initial  grain  size  distribution  affects  the  evolution 
behavior  within  the  thermokinetic  model.  In  applications  of  the  thermokinetic  model,  the  scalar 
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value  for  initial  grain  size  represents  the  average  grain  size.  However,  it  is  well  known  that  most 
polycrystalline  materials  have  a  lognormal  distribution  of  grain  sizes.  To  model  the  effect  of  a 
normal  and  lognormal  grain  size  distribution,  a  Monte  Carlo  analysis  is  applied  to  the 
thermokinetic  model,  whereby  the  initial  grain  size  was  randomly  chosen  1,000  times  from  these 
2  distributions  (with  a  2  nm  standard  deviation).  Figure  19  shows  the  results  of  these  calculations 
for  a)  a  normal  distribution  and  b)  a  lognormal  distribution  with  the  modified 
temperature-dependent  thermokinetic  model.  In  these  plots,  the  shaded  regions  denote  the  span  of 
the  grain  growth  curves  from  these  1,000  samples  for  each  temperature  (700,  800,  900,  and 
1,000  °C).  The  lines  within  the  distributions  denote  the  median  behavior  (thick  line),  the  10th 
percentiles  of  the  distribution  (dashed  lines),  and  the  minimum/maximum  values  of  the 
distribution  (bounding  lines  of  shaded  region).  Interestingly,  at  short  times,  the  initial  grain  size 
distribution  rapidly  approaches  the  average  grain  size  behavior  dave-  For  instance,  for  smaller 
grain  sizes  {d  <  dave),  there  is  a  larger  driving  force  for  grain  growth.  As  a  result,  these  grain  sizes 
rapidly  converge  to  the  mean  behavior.  For  larger  grain  sizes  (d  >  dave),  there  is  a  smaller  driving 
force  for  grain  growth  and  so  these  grain  sizes  very  slowly  evolve  towards  the  stabilized  grain 
size.  This  behavior  can  also  be  inferred  from  the  evolution  of  the  standard  deviation  of  these 
distributions  (Fig.  20).  Notice  that  the  standard  deviation  for  all  temperatures  and  both 
distributions  starts  at  2  nm  and  rapidly  decreases  until  it  reaches  a  linear  regime  (on  a  logarithmic 
scale,  though).  So,  with  respect  to  Figs.  19  and  20,  while  the  shaded  region  in  Fig.  19  reflects  a 
broad  initial  grain  size  distribution,  the  majority  of  the  distribution  rapidly  converges  to  the  mean 
behavior  while  grain  sizes  that  are  far  greater  than  dave  slowly  converge.  At  long  times,  the  initial 
distribution  of  grain  sizes  will  approach  a  single  stabilized  grain  size.  These  observations  are  true 
for  both  of  the  distributions. 
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(a) 


(b) 


Fig.  19  Grain  size  distribution  evolution  as  a  function  of 

annealing  time  for  a)  a  normal  distribution  and  b)  a 
lognormal  distribution.  The  lines  within  the 
distributions  denote  the  median  behavior  (thick  line), 
the  10th  percentiles  of  the  distribution  (dashed  lines), 
and  the  minimum/maximum  values  of  the 
distribution  (bounding  lines  of  shaded  region). 
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Standard  Deviation  in  Grain  Size  (nm) 


Fig.  20  Evolution  of  the  standard  deviation  of  the  grain  size  distribution  as  a 
function  of  annealing  time 
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5.  Discussion 


This  study  examined  the  thermodynamie  and  kinetie  parameters  of  the  as-formulated 
thermokinetie  model  and  proposed  extensions  to  the  formulation  to  eapture  the  temperature 
dependenee  of  grain  growth  within  this  model  formulation.  The  present  predictive  model  captures 
the  influence  of  temperature  better  than  previously-implemented  versions  of  the  thermokinetie 
model,  which  could  only  be  fit  to  existing  data  at  a  particular  temperature  and  had  little  predictive 
capability  for  grain  growth  outside  of  the  initial  conditions  (temperature,  concentration,  etc.).  For 
example,  compare  the  predictive  capability  of  the  original  model  in  Fig.  14  with  the  predictive 
capability  of  the  extended  model  in  Figs.  17  and  18.  As  further  illustration  of  this  difference. 

Fig.  21  shows  the  thermokinetie  model  a)  with  and  b)  without  the  temperature  dependence  in  the 
thermodynamic/kinetic  parameters.  In  the  case  of  the  original  model  (Fig.  21b),  the  temperature 
could  only  be  extended  ±50  °C  from  the  temperature  used  to  fit  the  model;  the  difference  in  the 
low  and  high  temperatures  (750  and  850  °C,  respectively)  with  this  model  resulted  in  a  difference 
in  grain  size  of  approximately  400%  after  1  h.  In  the  case  of  the  extended  model  herein  (Fig.  21a), 
the  same  temperature  increase  (AT  =100  °C)  shows  an  increase  of  approximately  10%  in  grain 
size  after  the  same  time  interval  (1  h),  in  agreement  with  experimental  results.  The 
experimentally-measured  XRD  data  over  a  400  °C  temperature  range  is  shown  as  black  diamonds 
in  Fig.  21  to  illustrate  the  magnitude  of  the  variation  in  grain  size  with  a  temperature  difference 
that  is  4  times  that  plotted  for  the  models  in  Fig.  21. 

Future  extensions  could  look  at  extending  this  thermokinetie  model  to  include  concentration 
dependence  as  well.  In  this  respect,  not  only  would  concentration  potentially  affect  the 
thermodynamic  and  kinetic  parameters,  but  a  number  of  material  parameters  would  also  be 
similarly  affected.  Additionally,  concentration  dependence  may  cause  problems  in  regimes  where 
rampant  grain  growth  exists  (Fe-0%  Zr)  or  intermetallic  formation  may  become  an  issue 
(Fe-10%  Zr).  Furthermore,  there  is  a  lot  of  inherent  noise  (and  uncertainty)  in  the  XRD 
measurements.  Nondestructive  techniques  that  can  dynamically  capture  the  evolution  of  grain  size 
with  time  with  less  uncertainty  would  certainly  improve  the  predictive  capability  of  the 
thermokinetie  grain  growth  model.  The  ability  to  predict  grain  growth  in  nanocrystalline  materials 
based  on  both  thermodynamic  and  kinetic  approaches  may  better  enable  high  temperature 
processing  and  consolidation  of  these  materials. 
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(a)  (b) 


Fig.  21  Evolution  of  the  average  grain  size  as  a  function  of  annealing  time  for  nanocrystalline  Fe-4%  Zr 
alloy  annealed  at  various  temperatures  in  the  range  750-850  °C.  Experimentally-measured  XRD 
data  for  a  400  °C  range  is  shown  (black  diamonds)  along  with  the  thermokinetic  model  a)  with 
and  b)  without  the  temperature  dependence. 


6.  Conclusions 


In  this  work,  a  thermokinetic  model  for  grain  growth  was  applied  to  experimental  XRD 
measurements  from  nanocrystalline  Fe-Zr  in  an  effort  1)  to  understand  the  influence  of 
thermodynamic,  kinetic,  and  material  parameters  in  the  model;  and  2)  to  extend  the  thermokinetic 
model  by  incorporating  temperature  dependence.  The  model  performs  well  for  the  grain  boundary 
saturated  case  in  the  nanocrystalline  Fe-Zr  system,  where  Zr  segregate  to  the  Fe  GBs  and 
thermodynamically/kinetically  reduce  the  driving  force  for  grain  growth.  In  this  work,  a 
sensitivity  analysis  of  parameters  (Monte  Carlo  global  sensitivity  analysis  with  10,000 
instantiations)  shows  that  only  2  independent  thermodynamic/kinetic  parameters  affect  the  grain 
size  evolution:  the  activation  energy  {Ql)  and  the  heat  of  segregation/grain  boundary  energy /grain 
boundary  saturated  excess  {HseghJ^bo,  respectively),  whose  effects  were  all  highly  correlated. 
Moreover,  the  predictive  capability  of  the  model  was  examined.  After  fitting  to  grain  growth  at  1 
temperature/concentration  condition,  the  model  performed  poorly  at  temperatures  as  close  as 
±50  °C  and  concentrations  as  close  as  ±0.5  at%.  This  model  was  then  extended  to  include  the 
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change  in  2  independent  thermodynamic/kinetic  parameters  as  a  function  of  temperature  to  model 
grain  growth  over  a  variety  of  temperatures  for  the  Fe-Zr  system,  in  agreement  with  the 
experimental  data  for  the  Fe-Zr  system.  Last,  a  Monte  Carlo  analysis  explored  the  effect  of  grain 
size  distribution  and  found  that  the  grain  size  distribution  rapidly  converged  to  the  mean  behavior. 
At  long  times,  the  grain  size  converged  to  a  single  stabilized  grain  size.  The  significance  of  this 
research  is  that  the  thermodynamic  and  kinetic  contributions  may  be  necessary  to  help  explain 
grain  growth  in  nanocrystalline  materials.  Furthermore,  this  model  formulation  can  be  extended 
to  understanding  how  grain  size  evolves  as  a  function  of  temperature  for  other  nanocrystalline 
systems. 
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Appendix  A.  Thermokinetic  Model  Validation  MATLAB  Script 


This  appendix  appears  in  its  original  form,  without  editorial  change. 
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%%  Thermokinetic  Grain  Growth  Model 


%  Mark  Tschopp,  U.S.  Army  Research  Laboratory,  2013-2014 
%  Email:  mark.a.tschopp@us.army.mil 

o, 

0 

%  The  thermokinetic  grain  growth  model  is  based  on  Reference  [1]  by  Chen 
%  and  colleagues.  This  is  the  basic  implementation  of  the  model  from  this 
%  publication.  A  basic  root  finding  method  is  used  to  calculate  the  grain 
%  size  as  a  function  of  time  for  this  analytical  model. 

o, 

0 

%  Reference: 

%  [1]  Chen,  Liu,  Yang,  Shen,  Acta  Materialia,  60  (2012)  4833. 

%  Clear  everything 
clear  all,  close  all 


%  Initialize  constants 
R  =  8.3145; 

V_m  =  7.1/ (100^3) ; 

M_0  =  0.8; 


%  J/mol  K;  ideal  gas  constant 
%  m3  mol-1;  Molar  volume 
%  m  mol  s-1  J-1;  Mobility  [1] 


%  Fitted  values: 
sigma_0  =  0.73; 

H_seg  =  89.3*1000; 
gamma_b0  =  1.25*10''(-5); 

Q_L  =  345.1*1000; 


%  J  m-2;  Grain  boundary  energy 
%  J  mol-1;  Heat  of  segregation 
%  mol  m-2;  Saturated  solute  excess 
%  J  mol-1;  Activation  energy  (lattice  diffusion) 


%  State  variables 
T  =  1073; 

x_0  =  0.04; 

D_0  =  5*10''  (-9)  ; 


%  K;  Temperature 
%  Zr  Concentration 
%  m;  initial  grain  size 


%  Assumptions 
g  amma_b  =  1 *  g amma_b  0 ; 
D_D0  =  1; 

D_L0  =  1; 

D  =  D_0; 


%  =gamma_b0  for  saturated  case;  solute  excess 
%  interface  diffusion 

%  lattice  diffusion  (ratio  gets  used) 


%%  Grain  Growth  Thermodynamics 

%  These  are  equations  13  and  14  [1]  for  the  saturated  grain  boundary  case 


sigma_l  =  sigma_0  -  gamma_b0* (R*T* (log (x_0 ) ) tH_seg) ; 
sigma_2  =  3*gamma_b0 ^2*R*T*V_m/x_0 ; 
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sigma_b  =  sigma_l  +  sigina_2*  (l/D)  ; 

disp  (  [  '  Sigma_l :  '  num2str  (  sigma_l )  '  J/m''2']) 

disp ( [ ' Sigma_2 :  '  num2str ( sigma_2 )  '  J/m']) 

disp ( [ ' Sigma_b :  '  num2str ( sigma_b)  '  J/m^2']) 

%%  Grian  Growth  Kinetics 

%  These  are  equations  20  and  21  [1]  for  the  saturated  grain  boundary  case 

Q_1  =  Q_L  +  R*T*log (D_D0/D_L0) -sigina_0/gamma_b0  t  (R*T* (log (x_0) ) tH_seg) ; 

Q_2  =  3*R*T*V_m*gamma_b0/x_0; 

Q_D  =  Q_l-Q_2/D; 

disp(['Q_l:  '  num2str (Q_l )  '  J/mol']) 

disp(['Q_2:  '  num2str (Q_2 )  '  J*m/mol']) 

disp(['Q_D:  '  num2str (Q_D)  '  J/mol']) 

%%  Stabilization  of  Nanocrystalline  alloys 

%  Stabilization  of  grain  growth  will  occur  when  GBs  are  saturated  and  there 
%  is  a  negligible  difference  between  $Q_L$  and  $Q_D$.  This  calculation  for 
%  D_star  shows  that  the  stabilized  grain  size  can  be  calculated  even  before 
%  the  "path"  to  the  stability  is  calculated.  Note  that  negative  D_star 
%  will  cause  problems  with  the  model. 

D_star  =  -sigma_2/sigma_l; 

disp ( [ ' D_star :  '  num2str (D_star*10^9)  '  nm']) 

%%  Implementation  of  Thermokinetic  Grain  Growth  Model 

%  This  section  implements  the  thermokinetic  grain  growth  model  (equation  24 
%  from  [1])  by  using  a  bisection  root  finding  method  to  solve  for  the  grain 
%  size,  d,  at  time  t. 

%  Insert  the  experimental  data 

t_exp  =  [60  15  10  7.5  2.5  0]  *  60;  %  s;  annealing  time 

D_exp  =  [26.5  20.5  21.3  19.5  14.3  7.4];  %  m;  grain  size  (experimental) 

%  To  match  Figure  3  in  [1],  the  end  time  is  set  to  1.2  hr 

t_end  =  1.2*60*60;  %  End  time  (in  seconds) 

t_a  =  [ ] ;  D_t  =  [ ] ; 

count  =  1; 

t_a (count)  =  0;  D_t (count)  =  D_0*10^9;  %  Initialize  time  and  grain  size 
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%  Main  loop  that  uses  numerical  root  finding  technique  to  solve  the 
%  analytical  thermokinetic  expression 
for  t  =  10 :t_end/100 :t_end  %  s 

count  =  count  +  1; 

xTol  =  O.OOle-10; 

xLower  =  D_0 ;  xUpper  =  D_star; 

LHS  =  0 . 5* ( (xLower+sigma_2 / sigma_l ) ^2- (D_0+sigma_2/ sigma_l ) ^2 ) - . . . 

(Q_2/ (R*T) +2*sigma_2/ sigma_l) * (xLower-D_0) + . . . 

(sigma_2^2/sigma_l''2  +  sigma_2/ sigma_l*Q_2/  (R*T)  +1/2*  (Q_2/  {R*T)  )  ^2)  *  .  .  . 
log { {sigma_l*xLower+sigma_2) / ( sigma_l *D_0+sigma_2 ) ) ;  % 

RHS  =  sigma_l*M_0*exp (-Q_l/ (R*T) ) *t; 
diff_S  =  LHS-RHS; 
y Lower  =  diff_S; 

xMid  =  (xLower  +  xUpper)/2.0; 

LHS  =  0 . 5*  (  (xMid+sigma_2  /  sigma_l )  ''2-  {D_0  +  sigma_2  /  sigma_l )  ''2)  -  .  .  . 

(Q_2 / (R*T) +2*sigma_2/ sigma_l ) * (xMid-D_0 ) + . . . 

(sigma_2^2/sigma_l''2  +  sigma_2/ sigma_l*Q_2/  (R*T)  +1/2*  (Q,_2/  {R*T)  )  ’^2)  *  .  .  . 
log{ (sigma_l*xMid+sigma_2) / (sigma_l*D_0+sigma_2) ) ;  % 

RHS  =  sigma_l*M_0*exp (-Q_l/ (R*T) ) *t; 
diff_S  =  LHS-RHS; 
yMid  =  diff_S; 

%  Calculate  left-hand-side  (LHS)  and  right-hand-side  (RHS)  of  equation 
%  and  iteratively  hone  in  on  the  grain  size  solution  (xMid) . 
while  (xUpper  -  xLower) /2.0  >  xTol 
if  yLower  *  yMid  >  0.0 
xLower  =  xMid; 

else 

xUpper  =  xMid; 

end 

xMid  =  (xLower  +  xUpper) /2.0; 

LHS  =  0 . 5*  (  (xMid+sigma_2  /  sigma_l )  ''2  -  (D_0  +  sigma_2/  sigma_l )  ''2)  -  .  .  . 

(Q_2/ (R*T) +2*sigma_2/ sigma_l) * (xMid-D_0) + . . . 

(sigma_2^2/ sigma_1^2+sigma_2/sigma_l*Q_2/ (R*T) +1/2* (Q_2/ (R*T) ) ^2) * . . . 
log ( ( sigma_l *xMid+sigma_2 ) / ( sigma_l *D_0+sigma_2 ) ) ;  %  Calculate  LHS 

RHS  =  sigma_l*M_0*exp (-Q_l/ (R*T))*t;  %  Calculate  RHS 

diff_S  =  LHS-RHS; 
yMid  =  diff_S; 

end 

if  xMid  >  D_star,  xMid  =  D_star;  end 
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t_a { count ) 


t; 


D_t  (count)  =  xMid*10''9; 


end 

%%  Plot  the  results 

%  Plot  the  experimental  data 

plot  (t_exp/3600 ,  D_exp,  '  db  ' ,  '  MarkerFaceColor  ' ,  '  b  '  ) /■  hold  on 

xlim ( [0  1.2]) 

%  Plot  the  thermokinetic  model  data 
plot (t_a/ 3600, D_t , ' -b ' , ' LineWidth ' , 2 ) 

axis  square 

legend { ' Fe-Zr  Expt ' , ' Thermokinetic  Model ' , ' Location ',[0.71  0.2  0.01  0.01]) 
xlabel { 'Annealing  Time  (hr)  '  ,  .  .  . 

'FontSize ' , 30,  ' FontName ',' Times '  ,  ' FontWeight '  ,  ' b ' )  ; 

ylabel ( ' Grain  Size  (nm)  '  ,  .  .  . 

'FontSize ', 30,  ' FontName ',' Times ' ,  ' FontWeight ',' b ') ; 

set (gca,  ' FontSize ', 24,  ' FontName ',' Times ' ,  ' LineWidth ', 2 ) ; 

set (get (gca,  ' Children ' ) ,  'Marker Size ' , 10) ; 
axis  (  [0  1.2  5  30] ) 

xlabel (get (get (gca, ' xlabel '),' String ') ,  ' FontSize ' , 30 ,  ' FontName ' , ' Times ' , 

' FontWeight ' , ' b ' ) ; 

ylabel (get (get (gca,  ' ylabel '),' String ') ,  ' FontSize  '  , 30 ,  ' FontName ' ,  ' Times ' , 

' FontWeight ' , ' b ' ) ; 

set (gca,  ' FontSize ', 24,  ' FontName ',' Times ' ,  ' LineWidth ', 2 ) ; 

title ( ' remove ' ) 
legend( 'boxoff ' ) 
axis  square 

set (gcf, ' Position ',[ 1  1  1200  1200]) 
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Intentionally  lelt  blank. 
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Appendix  B.  Thermokinetic  Model  Optimization  MATLAB  Script 


This  appendix  appears  in  its  original  form,  without  editorial  change. 
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%%  Thermokinetic  Grain  Growth  Model 

g, 

0 

%  Mark  Tschopp,  U.S.  Army  Research  Laboratory,  2013-2014 
%  Email:  mark.a.tschopp@us.army.mil 

o, 

0 

%  The  thermokinetic  grain  growth  model  is  based  on  Reference  [1]  by  Chen 
%  and  colleagues. 

o, 

o 

%  Reference: 

%  [1]  Chen,  Liu,  Yang,  Shen,  Acta  Materialia,  60  (2012)  4833. 

g, 

0 

%  In  this  code,  the  Fe-4%Zr  data  is  used  and  the  optimization  routine  is 
%  outlined  for  the  modified  thermokinetic  model  with  temperature  dependence 
%  of  thermodynamic/kinetic  parameters.  This  is  performed  in  a  series  of 
%  optimization  routines,  to  ensure  convergence  of  the  model. 

%  Clear  everything 
clear  all,  close  all 

%  Import  data  from  Excel  file 

filename  =  'Fe  Zr  hot  XRD  complete  data  set.xls'; 

[num, txt , raw]  =  xlsread (filename) ; 

%  'inc'  specifies  which  dataset  to  use,  i.e.,  in  this  case  it  is  the  4%Zr 

%  and  700,  800,  900,  and  1000  degrees  Celsius 

inc  =  [0000,  0000,  0000,  1111,  0000]; 

temp  =  [400  500  600  700,  700  800  900  1000,  700  800  900  1000,  .  .  . 

700  800  900  1000,  700  800  900  1000]+273; 
cone  =[0000,  0.005  0.005  0.005  0.005,  0.01  0.01  0.01  0.01,  ... 

0.04  0.04  0.04  0.04,  0.1  0.1  0.1  0.1]; 

%  Material  parameters 

R  =  8.3145;  %  J/mol  K;  ideal  gas  constant 

V_m  =  7.1/ (100^3);  %  m3  mol-1;  Molar  volume 

%  Use  these  two  independent  thermodynamic/thermokinetic  parameters 
hseg  =  zeros  (size (inc) ) ; 
ql  =  zeros (size (inc) ) ; 

%%  Optimization  I 

%  In  this  optimization,  we  are  solving  for  the  optimal  sigma_0  and  Q_L 
icount  =  0; 
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f  =  0; 

for  i  =  1:2:39 

imod  =  ( i+1 ) 12} 
if  inc (imod) 

icount  =  icount  +  1; 
t_expl ( : , icount )  =  num(:,i); 

D_expl  (  : , icount )  =  num ( : , i  +  1 )  '  ; 

Temp (icount)  =  temp (imod); 

Cone (icount)  =  cone (imod); 

T  =  temp (imod);  x_0  =  cone (imod); 
t_exp  =  num ( : , i ) ;  D_exp  =  num ( : , i+1 ) ' ; 

%  Select  starting  values 


sigma_0  =  0.73; 

"6 

J  m-2 ; 

H_seg  =  100*1000; 

g, 

o 

J  mol-1; 

Heat  of  segregation 

gamma_b0  =  1.25*10''(-5); 

"6 

mol  m-2; 

saturated  solute  excess 

Q_L  =  345.1*1000; 

o, 

o 

J  mol-1; 

lattice  diffusion 

o, 

o 

activation  energy 

save ( ' ob jfun_mod_exp_sigmaO .mat ' , ' t_exp ' , ' D_exp ' , ' T ' , ' x_0 ' , . . . 

'  H_seg '  ,  '  gainma_b0  '  )  ; 
xO  =  [sigma_0  Q_L]  ; 
options  =  optimset (' Display of f ')  ; 
xsol  =  fminsearch ( @obj fun_mod_sigmaO , xO , options ) ; 

sigma_0  =  xsol(l);  %  J  mol-1;  Heat  of  segregation 

Q_L  =  xsol (2);  %  J  mol-1;  lattice  diffusion  activation  energy 

disp  (  [ ' xsol ( 1 )  :  '  num2str (xsol ( 1 ) ) ] ) 

disp ( [ ' xsol (2 ) :  '  num2str (xsol (2) ) ] ) 

f  =  f  +  ob jfun_mod_sigmaO (xsol)  ; 

sigmaO (imod)  =  xsol(l); 
ql(imod)  =  xsol (2); 

%  Are  initial  parameters  giving  a  positive  D_star? 
sigma_l  =  sigma_0  -  gamma_b0* (R*T* (log (x_0) ) +H_seg) ; 
sigma_2  =  3*gamma_b0^2*R*T*V_m/x_0; 

D_star  =  -sigma_2/sigma_l ; 

if  D_star  <  0,  stop;  end 

end 

end 
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%  'f  is  a  measure  of  the  goodness  of  the  model  fit,  e.g.,  a  root  mean 
%  square  error  or  the  like 

disp ( [ ' f :  '  num2str (ob jfun_mod_sigmaO (xsol) ) ] ) 

%%  Optimization  II 

%  Now  use  sigma_0  to  get  good  estimates  of  what  the  constants  describing 
%  the  temperature  dependence  are.  In  this  case,  Q_L  is  constant. 

n  =  find (logical (inc) ) ; 

sigma_0_cl  =  mean { { sigma 0 (n { 2 : end) )  -  sigmaO (n ( 1 ) ) ) . /temp (n { 2 : end) ) ) ; 
sigma_0_c0  =  sigmaO(n(l))  -  sigma_0_cl*temp (n (1) ) ; 

Q_L  =  mean (ql  (n)  ) ; 

save ( ' ob j f un_mod2_exp_sigma0 . mat ' ,  ' t_expl ' ,  ' D_expl ' ,  ' Temp ' ,  ' Cone ' ,  . .  . 

' H_seg ' , ' gamma_bO ' ,  ' Q_L ' ) ; 

xO  =  [sigma_0_c0  sigma_0_cl ] ; 
options  =  optimset {' Display of f ')  ; 
xsol  =  fminsearch(@objfun_mod2_sigma0,x0,options); 

disp  ([' xsol ( 1 )  :  '  num2str (xsol ( 1 ) ) ] ) 

disp ( [ ' xsol (2 ) :  '  num2str (xsol (2 ) ) ] ) 

sigma_0_c0  =  xsol(l); 
sigma_0_cl  =  xsol(2); 

disp ( [ ' f :  '  num2str (ob jfun_mod2_sigma0 (xsol) ) ] ) 

%%  Optimization  III 

%  Now  use  Q_L  to  get  good  estimates  of  what  the  constants  describing 
%  the  temperature  dependence  are.  In  this  case,  sigma_0  is  constant.  In 
%  this  routine,  multiple  starting  positions  are  utilized  to  find  the 
%  minimum  error  solution  to  the  thermokinetic  equations. 

Q_L_cO  =  Q_L;  Q_L_cl  =  1; 
best  =  100; 

for  i  =  0.7:0. 1:1. 3 
for  j  =  0  :  0 . 5  :  5 

Q_L_c  0  =  i *  Q_L ; 

Q_L_cl  =  10"j; 
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save (' objfun_mod3_exp_sigmaO_Arrhenius .mat ' , 't_expl', 'D_expl', ... 

' Temp ' , ' Cone ' , ' H_seg ' , ' gamma_bO ' , ' sigma_0_c0 ' , ' sigma_0_cl ' ) ; 

X  0  =  [ Q_L_c  0  Q_L_c 1 ] ; 

options  =  optimset (' Display of f ')  ; 

xsol  =  fminsearch ( @obj f un_mod3_sigmaO_Arrhenius , xO , options )  ; 

disp ( [ ' xsol  ( 1 )  ;  '  num2str (xsol { 1 ) ) ] ) 

disp ( [ ' xsol (2 ) :  '  num2str (xsol (2) ) ] ) 

Q_L_cO  =  xsol(l); 

Q_L_cl  =  xsol (2); 

disp(['f:  '  num2str (ob jfun_mod3_sigmaO_Arrhenius (xsol) ) ] ) 
f  =  ob jfun_mod3_sigmaO_Arrhenius (xsol)  ; 
if  f  <=  best 

a  =  0  L  cO ; 
b  =  Q_L_c 1 ; 
best  =  f; 

end 

end 

end 

0  L  cO  =  a; 

Q_L_cl  =  b; 

%%  Optimization  IV 

%  In  this  optimization,  the  previous  solutions  for  the  temperature 
%  dependent  constants  are  used  as  starting  points  for  the  final 
%  optimization.  Without  these  different  steps,  the  optimization  may  not 
%  converge  to  a  valid  solution. 

save ( ' ob jfun_mod4_exp_sigmaO_Arrhenius .mat ' ,  ' t_expl ' ,  ' D_expl ' ,  ' Temp '  ,  .  .  . 

' Cone ' , ' H_seg ' , ' gamma_bO ' ) ; 
xO  =  [sigma_0_c0  sigma_0_cl  Q_L_cO  Q_L_cl]; 

options  =  optimset (' Display of f ')  ; 

xsol  =  fminsearch ( @obj f un_mod4_sigmaO_Arrhenius , xO , options ) ; 

disp ([' xsol ( 1 ) :  '  num2str (xsol ( 1 ) ) ] ) 

disp  (  [ ' xsol (2 )  :  '  num2str (xsol (2 ) ) ] ) 

disp ( [ ' xsol ( 3 ) :  '  num2str (xsol (3) ) ] ) 

disp ( [ ' xsol ( 4 ) :  '  num2str (xsol (4) ) ] ) 
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sigma_0_c0  =  xsol{l); 
sigma_0_cl  =  xsol{2); 

Q_L_cO  =  xsol(3); 

Q_L_cl  =  xsol(4); 

disp(['f:  '  num2str (ob jfun_mod4_sigmaO_Arrhenius (xsol) ) ] ) 

%  It  is  apparent  from  the  value  of  'f  (goodness  of  fit)  that  the  fit  is 
%  not  as  good  as  individually  fitting  the  model  to  each  specific 
%  temperature/concentration  combination.  But,  the  following  fit  better 
%  depicts  the  temperature  dependence  aspects  than  any  of  these  individual 
%  fits. 


%%  Plot 

%  Plot  the  experimental  results  with  the  model  results 

xO  =  [sigma_0_c0  sigma_0_cl  Q_L_cO  Q_L_cl]; 
hold  on 


%  Plot  the  thermokinetic  model  curves 
for  i  =  1 : length (Temp) 

[t_a,  D_t ]  =  ob j f un_mod4_plot_sigmaO_Arrhenius (xO , i ) ; 


if  i  ==  1 

hi  =  plot  (  [  0 
elseif  i  ==  2 

h2  =  plot  (  [  0 
elseif  i  ==  3 

h3  =  plot  (  [  0 
elseif  i  ==  4 

h4  =  plot  (  [  0 

end 

end 


t_a] ,  [D_expl(l,i)  D_t 
t_a] ,  [D_expl(l,i)  D_t 
t_a] ,  [D_expl(l,i)  D_t 
t_a] ,  [D_expl(l,i)  D_t 


]  ,  ' -r ' ,  ' LineWidth ' , 3 ) ; 
]  ,  ' -g ' ,  ' LineWidth ' , 3 ) ; 
] , ' -c ' , ' LineWidth ' , 3 ) ; 
]  ,  ' -b ' ,  ' LineWidth ' , 3 ) ; 


%  Plot  the  experimental  data 
count  =  0; 
for  i  =  1:2:39 

X  =  num ( : , i ) ; 
y  =  num ( : , i+1 ) ; 
imod  =  (i  +  1 )  12} 
if  inc (imod) 

count  =  count  +  1; 
if  count  ==  1 

hi  =  plot(x,y,  'or',  ' Markers ize ' , 8,  ' MarkerFaceColor ' ,  ' r ' ) ; 
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elseif  count  ==  2 

h2  =  plot(x,y,  'og',' Markers ize '  ,  8 ,  ' MarkerFaceColor '  ,  '  g  '  )  ; 
elseif  count  ==  3 

h3  =  plot (x, y,  ' sc MarkerSize 8 MarkerFaceColor c ')  ; 
elseif  count  ==  4 

h4  =  plot(x,y,  'sb',' MarkerSize ' , 8 ,  ' MarkerFaceColor ' ,  ' b ' )  ; 

end 

if  count  >=  1 

ylim( [5  35] ) 
xlim( [-100  6000] ) 
if  1  ==  7 

ylim([0  600]) 

legend] [hi  h2  h3  h4],'nc  Fe  400\circ { C }  ' ,  ' nc  Fe  500\circ{C} '  ,  . 
' nc  Fe  600\circ [ C } ' ,  ' nc  Fe  7  00\circ [ C  }  '  ,  .  .  . 

'Location' , [0.76  0.81  0.01  0.01]); 
elseif  i  ==  15 

legend([hl  h2  h3  h4],'nc  Fe-0.5%Zr  7 0 0 \circ [ C } ' , . . . 

' nc  Fe-0.5%Zr  800\circ { C } ' , ' nc  Fe-0.5%Zr  90 0 \circ [ C } ' , . . . 

' nc  Fe-0.5%Zr  lOOOXcirc (C } ' ,  .  .  . 

'Location' ,  [0 . 73  0.21  0.01  0.01]); 
elseif  1  ==  23 

legend] [hi  h2  h3  h4],'nc  Fe-l%Zr  700\circ[C} ' , . . . 

'nc  Fe-l%Zr  800\circ { C } ' , ' nc  Fe-l%Zr  900\circ{C} ' , . . . 

'nc  Fe-l%Zr  1000\circ{C}  '  ,  .  .  . 

'Location' , [0.73  0.21  0.01  0.01]); 
elseif  1  ==  31 

legend ][hl  h2  h3  h4],'nc  Fe-4%Zr  700\circ[C}',... 

' nc  Fe-4%Zr  800\circ ( C } ' ,  ' nc  Fe-4%Zr  900\circ{C}  '  ,  .  .  . 

' nc  Fe-4%Zr  1000\circ{C}  ',  .  .  . 

'Location' ,  [0 . 73  0.21  0.01  0.01]); 
elseif  1  ==  39 

legend]  [hi  h2  h3  h4],'nc  Fe-10%Zr  70 0 \circ [ C } ' ,  .  .  . 

'nc  Fe-10%Zr  800\circ { C } ' , ' nc  Fe-10%Zr  900\circ [ C } ' , . . . 

'nc  Fe-10%Zr  1000\circ{C}  '  ,  .  .  . 

'Location' ,  [0.73  0.21  0.01  0.01]); 

end 

legend ] ' boxof f ' ) 

X label ] ' Time  ]  s )  ' ,  .  . . 

'FontSize ' , 30,  ' FontName ',' Times ' ,  ' FontWeight ' , ' b ' ) ; 

y label ['Grain  Size  ]nm) ' , . . . 

'FontSize ',  30,  ' FontName ',' Times ' ,  ' FontWeight ',' b ') ; 

set ]gca,  'FontSize ', 24,  ' FontName ',' Times ' ,  ' LineWidth ' , 2 , . . . 

' TickDir '  ,  ' out ' ,  ' XMinorTick ' ,  ' on ' ,  ' YMinorTick ' ,  ' on ' ) ; 
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set (get (gca,  'Children ')  ,  'Marker Size ' , 8,  ' LineWidth ' , 2,  .  .  . 

' MarkerEdgeColor '  ,  'k') 
axis  square 

xlabel (get (get (gca, ' xlabel '),' String ') ,  ' FontSize ' , 30,  ... 

' FontName ' , ' Times ' ,  ' FontWeight ' , 'b'); 

ylabel (get (get (gca,  ' ylabel '),' String ') ,  'FontSize ', 30,  ... 

' FontName ' ,  ' Times ' ,  ' FontWeight ' ,  ' b ' )  ; 

set (gca,  'FontSize ',  24,  ' FontName ',' Times ' ,  ' LineWidth ', 1 ) ; 

title ( ' remove ' ) 
legend ( ' boxof f ' ) 
axis  square 

set (gcf,  ' Position ',[  1  1  1200  1200]) 

end 

end 

end 
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Appendix  C.  Thermokinetic  Model  Optimization  Function  MATLAB  Script 


This  appendix  appears  in  its  original  form,  without  editorial  change. 
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%%  Thermokinetic  Grain  Growth  Model 


%  Mark  Tschopp,  U.S.  Army  Research  Laboratory,  2013-2014 
%  Email:  mark.a.tschopp@us.army.mil 

o, 

0 

%  The  thermokinetic  grain  growth  model  is  based  on  Reference  [1]  by  Chen 
%  and  colleagues. 

o, 

o 

%  Reference: 

%  [1]  Chen,  Liu,  Yang,  Shen,  Acta  Materialia,  60  (2012)  4833. 

o, 

0 

%  This  is  the  optimization  function  used  for  the  final  optimization: 

function  [t_a,  D_t ]  =  ob jfun_mod4_plot_sigmaO_Arrhenius (x, c) 

sigma_0_c0  =  x(l); 
sigma_0_cl  =  x(2); 

Q_L_c  0  =  X ( 3 ) ; 

Q_L_cl  =  x(4)  ; 


R  =  8.3145; 

g, 

0 

J/mol  K;  ideal  gas  constant 

V_m  =  7.1/ (100^3)  ; 

o_ 

o 

m3  mol-1;  Molar  volume 

o 

II 

o 

CO 

o_ 

0 

m  mol  s-1  J-1;  Interface  mobility 

D_D0  =  1; 

o_ 

0 

interface  diffusion 

D_L0  =  1; 

g_ 

0 

lattice  diffusion 

t_expl  =  [];  D_expl  = 

[] ; 

load ( ' ob jfun_mod4_exp_sigmaO_Arrhenius .mat ' )  ; 
i  =  c; 

T  =  Temp { i ) ; 
x_0  =  Cone ( i ) ; 
t_exp  =  t_expl ( : , i) ; 

D_exp  =  D_expl ( :  ,  i)  '  ; 

sigma_0  =  sigma_0_c0  t  sigma_0_cl*T; 

Q_L  =  Q_L_cO*exp  (-sqrt  (Q_L_cl''2)  /  (R*T)  )  ; 

sigma_l  =  sigma_0  -  gamma_bO* {R*T* (log {x_0 ) ) tH_seg) ; 
sigma_2  =  3*gamma_b0 ^2*R*T*V_m/x_0 ; 

Q_1  =  Q_L  +  R*T*log (D_D0/D_L0) -sigma_0/gamma_b0  t  (R*T* (log (x_0) ) tH_seg) ; 
Q_2  =  3*R*T*V_m*gamma_b0/x_0; 
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D_star 


=  -sigma_2/sigma_l; 

clear  t_a  D_t 
count  =  0; 

t_plot  =  t_exp; 

D_0  =  D_exp (1) *10" (-9)  ; 

for  ti  =  1 : length (t_plot ) 
t  =  t_plot (ti)  ; 

count  =  count  +  1; 
xTol  =  O.OOle-10; 
xLower  =  D_0 ; 
xUpper  =  D_star; 

LHS  =  0 . 5* ( (xLower+sigma_2 / sigina_l ) "2- (D_0+sigma_2 / sigma_l ) "2 ) - . . . 
(Q_2/ (R*T) +2*sigma_2/ sigma_l) * (xLower-D_0) + . . . 

(sigma_2"2/sigma_l"2+sigma_2/ sigma_l*Q_2/ (R*T) +1/2* (Q_2/ {R*T) ) "2) * 
log{ (sigina_l*xLower+sigma_2) / ( sigma_l *D_0+sigma_2 ) ) ;  % 

RHS  =  sigina_l*M_0*exp  (-Q_l/ (R*T)  )  *t; 
diff_S  =  LHS-RHS; 
y Lower  =  diff_S; 

xMid  =  (xLower  +  xUpper)/2.0; 

LHS  =  0 . 5* ( (xMid+sigma_2 / sigma_l ) "2- {D_0+sigma_2 / sigina_l ) "2 ) - . . . 
(Q_2 / (R*T) +2*sigma_2/ sigma_l ) * (xMid-D_0 ) + . . . 

(sigma_2"2/sigma_l"2+sigma_2/ sigma_l*Q_2/ (R*T) +1/2* (Q_2/ {R*T) ) "2) * 
log{ (sigina_l*xMid+sigma_2) / (sigma_l*D_0+sigma_2) ) ;  % 

RHS  =  sigma_l*M_0*exp (-Q_l/ (R*T) ) *t; 
diff_S  =  LHS-RHS; 
yMid  =  diff_S; 

iters  =  0; 

while  (xUpper  -  xLower) /2.0  >  xTol 
iters  =  iters  +  1; 
if  yLower  *  yMid  >  0.0 
xLower  =  xMid; 

else 

xUpper  =  xMid; 

end 

xMid  =  (xLower  +  xUpper) /2.0; 

LHS  =  0 . 5* ( (xMid+sigma_2 /sigma_l ) "2- (D_0+sigma_2/sigma_l ) "2 ) - . 
(Q_2/  (R*T)  +2*sigina_2/ sigma_l)  *  (xMid-D_0)  +  .  .  . 
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{sigina_2^2/ sigma_1^2+sigma_2/sigma_l*Q_2/ {R*T) +1/2* (Q_2/ (R*T) ) ^2) * . . . 
log ( ( sigma_l *xMid+sigma_2 ) / ( sigma_l *D_0+sigma_2 ) ) ;  % 

RHS  =  sigma_l*M_0*exp {-Q_l/ (R*T) ) *t; 
diff_S  =  LHS-RHS; 
yMid  =  diff_S; 

end 

If  xMid  >  D_star,  xMid  =  D_star;  end 
t_a{ count)  =  t; 

D_t  (count)  =  xMid*10''9; 


end 

end 
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